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Abstract 


From  1999  through  2001  the  University  of  Michigan  carried  out  an  experimental  study  of 
the  high  Reynolds  number  turbulent  flow  over  a  hydrofoil.  The  high  Reynolds  number  foil,  or 
Hifoil  project  was  carried  out  in  three  phases  at  the  U.S.  Navy’s  William  B.  Morgan  Large 
Cavitation  Channel  (LCC),  the  world’s  largest  variable  pressure  recirculating  water  tunnel. 
Several  standard  flow  measurements  were  made  over  the  course  of  this  project,  including  surface 
static  pressure  field  measurements.  Static  pressure  was  measured  through  the  use  of  surface 
static  pressure  taps  drilled  into  the  foil  model  at  various  chordwise  and  spanwise  locations. 
Pressure  data  was  collected  by  a  custom  designed  and  built  static  pressure  acquisition  system, 
which  could  automate  the  process  if  so  desired.  The  foil  model  was  a  modified  NACA  16  shape 
with  a  chord  of  7  feet  and  a  span  of  10  ft.  The  baseline  foil  was  tested  at  0,  -1  and  +1  degrees 
angle  of  attack  measure  relative  to  the  foil’s  flat  pressure  side.  Additional  data  was  collected  on 
the  foil  with  a  modified  trailing  edge  shape  again  at  0  degrees  angle  of  attack.  Data  was 
collected  at  flow  speeds  ranging  from  1.5  to  18.3  m/s  and  at  temperatures  ranging  from  22  to 
40  °C  (70  to  104  °F),  resulting  in  a  chord  based  Reynolds  number  range  of  3  to  60  million.  The 
Reynolds  number  range  extends  from  the  existing  model  scale  up  to  full  scale  for  some  naval 
applications.  Data  is  presented  here  in  the  form  of  normalized  nondimensional  pressure  as  a 
function  of  chordwise  location  represented  in  percentage  of  chord  at  all  the  various  geometric 
and  flow  speed  test  conditions.  In  addition  a  numeric  integration  of  the  pressure  curves  was 
performed  to  calculate  the  pressure  lift  generated  by  the  foil  at  the  various  test  conditions.  This 
data  is  also  presented  in  the  form  of  lift  coefficient  as  a  function  of  foil  configuration  and 
Reynolds  number.  The  pressure  data  presented  shows  the  expected  increase  in  lift  with 
increased  angle  of  attack  and  little  to  no  Reynolds  number  dependence. 
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Introduction 


From  the  fall  of  1999  through  the  summer  of  2001  an  experiment  was  conducted  by  the 
University  of  Michigan  to  investigate  the  turbulent  flow  over  a  large  hydrofoil  shape  relevant  to 
the  U.S  Navy.  The  experiments  were  performed  at  a  range  of  Reynolds  numbers  that  reached 
full  scale  for  some  actual  naval  applications.  All  experiments  associated  with  this  project  were 
conducted  at  the  U.S.  Navy’s  William  B.  Morgan  Large  Cavitation  Channel  (WBM-LCC),  the 
world’s  largest  variable  pressure  recirculating  water  tunnel.  The  High  Reynolds  Number  Foil,  or 
Hifoil  project  was  carried  out  for  a  variety  of  reasons.  One  important  goal  of  the  experiment  was 
to  collect,  under  controlled  experimental  settings,  a  high  Reynolds  number  hydrofoil  data  set  for 
the  U.S.  Navy  which  could  bridge  the  gap  between  the  existing  model  scale  data  and  the  full 
scale. 

The  data  set  produced  by  the  Hifoil  project  will  be  used  in  part  to  prove  that  what  is 
understood  about  this  type  of  a  flow  field  at  lower  Reynolds  numbers  actually  does  scale  up  to 
higher  Reynolds  numbers  in  a  predictable  manner.  As  a  result,  it  was  important  for  the  data 
collected  during  this  project  to  be  of  the  kind  appropriate  for  use  in  validating  scaling  laws  and 
computational  models  that  were  developed  based  on  the  existing  model  scale  data.  The  Navy 
obviously  knows  a  great  deal  about  the  flow  over  various  foil  shapes  on  the  model  scale, 
however  there  is  very  little  data  concerning  this  type  of  a  flow  over  foil  shapes  at  the  full  scale 
due  to  the  costs  and  complications  associated  with  collecting  high  Reynolds  number  data. 
Because  the  Hifoil  project  will  be  able  to  achieve  essentially  full  scale  chordwise  Reynolds 
numbers,  one  of  the  major  purposes  of  the  project  was  to  provide  information  that  can  be  readily 
used  by  the  Navy  in  Reynolds  number  scaling  studies. 
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To  this  end,  several  relatively  simple  measurements  were  made  during  the  Hifoil  project 
in  order  to  compare  the  flow  field  to  that  of  similar  experiments  conducted  at  lower  Reynolds 
numbers.  Among  these  were  trailing  edge  velocity  measurements,  foil  acceleration  levels  and 
both  dynamic  and  static  surface  pressure  measurements.  The  foil  surface  pressure  gradients 
measured  over  the  course  of  the  experiment  were  of  particular  importance.  Surface  pressure 
gradients  are  used  to  help  understand  various  aspects  of  a  flow  field,  such  as  the  amount  of  lift 
generated.  Static  pressure  is  a  simple  but  critical  measurement  that  must  be  made  in  any  foil 
experiment.  As  a  result  there  is  a  vast  amount  of  preexisting  static  pressure  data  from  various 
different  foil  shapes  at  various  Reynolds  numbers.  This  fact  made  the  static  pressure  data  one  of 
the  critical  data  sets  collected  during  the  Hifoil  project  because  it  offers  a  direct  comparison  to 
data  taken  during  virtually  every  other  foil  experiment  done  in  the  recent  past.  To  see  recent 
applications  of  static  pressure  data  see  Astolfi  [1],  Niewald  [7,8],  Wood  [10]  and  Zandieh  [11]. 

In  order  to  measure  the  static  pressure  during  the  Hifoil  project,  several  issues  had  to  be 
addressed.  Essentially  all  static  pressure  measurements  made  in  recent  memory  were  made  using 
small  surface  pressure  tap  holes  drilled  orthogonal  to  the  surface  of  the  experimental  model.  The 
same  method  was  implemented  for  the  Hifoil  project.  However  the  number,  size  and  location  of 
the  pressure  taps  to  be  used  had  to  be  determined  while  giving  careful  consideration  to  any  errors 
or  disturbances  that  the  taps  may  be  introducing  into  both  the  pressure  measurements  and  the 
flow  field  in  general.  Next,  the  sensors  to  be  used  had  to  be  selected  carefully  while  considering 
the  full  range  of  pressures  that  the  foil  was  likely  to  experience  over  the  course  of  the 
experiment.  Lastly,  a  pressure  acquisition  system  had  to  be  designed  and  implemented  which 
could  measure  the  fluid  pressure  in  the  various  taps  accurately  and  in  a  reasonable  amount  of 


time.  The  time  required  to  make  each  measurement  on  the  Hifoil  project  was  always  an 
important  factor  due  to  the  high  operating  costs  of  the  experimental  facility. 

A  great  deal  of  work  has  been  done  concerning  the  process  of  measuring  static  pressure 
in  a  flow  with  wall  pressure  taps.  The  work  of  Franklin  and  Wallace  [6],  Ducruet  [4,5]  and 
Zogg  [12],  among  others  has  provided  a  set  of  guidelines  to  be  observed  when  measuring  static 
pressure  with  wall  taps.  First,  the  taps  must  be  drilled  in  the  surface  of  the  model  orthogonal  to 
the  flow.  Second,  the  taps  must  have  as  clean  and  sharp  of  edges  as  possible.  The  tap  itself 
should  be  a  clean  and  flat  hole  in  the  surface  of  the  model  with  no  rounding  of  its  edges  or 
distortions  in  the  surface  of  the  model  upstream  or  downstream  of  the  tap.  Any  of  these 
problems  will  create  undesired  errors  in  the  pressure  measurements.  However,  all  the  pressure 
tap  work  done  has  shown  that  errors  in  the  pressure  measurements  are  unavoidable  because  the 
presence  of  the  pressure  tap  on  the  model  will  alter  the  flow  over  the  tap  slightly  no  matter  how 
well  it  has  been  machined  into  the  surface  of  the  model.  Benedict  compiled  data  from  various 
pressure  tap  studies  and  created  a  correlation  for  this  unavoidable  static  pressure  tap  hole  error. 

[3]  This  correlation  is  function  of  hole  diameter,  flow  speed,  wall  shear  stress,  fluid  viscosity 
and  fluid  density.  It  was  developed  using  the  data  from  four  different  pressure  tap  studies, 
including  Franklin  and  Wallace’s  data.  As  a  result  this  correlation  can  be  taken  to  be  fairly 
accurate,  assuming  the  static  pressure  taps  are  properly  machines  as  earlier  described.  This 
pressure  tap  theory  was  observed  in  the  designing  of  the  static  pressure  acquisition  system  for  the 
Hifoil  project,  in  the  selection  of  the  transducers  and  other  components  used  on  the  project,  in  the 
machining  of  the  pressure  taps  into  the  foil  and  in  the  calculation  of  the  uncertainty  associated 
with  the  static  pressure  measurements  made  over  the  course  of  the  project. 
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Throughout  the  Hifoil  project  surface  pressure  measurements  were  made  at  various 
chordwise  locations  on  both  the  suction  and  pressure  side  of  the  foil.  These  measurements  were 
repeated  at  various  speeds  and  for  various  angles  of  attack.  Based  on  these  pressure 
measurements  pressure  lift  was  then  calculated  at  all  the  various  conditions  in  order  to  get  a 
complete  picture  of  the  effect,  if  any,  that  high  Reynolds  number  has  on  lift.  Errors  arose  in  the 
measurements  due  to  various  sources  such  as  transducer  error,  calibration  offset  errors  and  static 
pressure  tap  hole  errors.  However,  all  these  sources  of  error  were  anticipated,  quantified  in  some 
reasonable  manner  and  then  either  reported  or  corrected  for  within  the  data. 


7 


Experimental  Procedure 


Project  Overview 

The  Hifoil  project  was  conducted  in  a  series  of  three  phases  at  the  WBM-LCC  over  the 
course  of  two  years.  A  series  of  experiments  were  conducted  on  a  modified  NACA  16  foil.  The 
modifications  made  to  the  base  NACA  16  foil  included  a  flat  pressure  side  to  simplify  the 
construction  of  the  model  and  an  increased  knuckle  on  the  trailing  edge  in  order  to  prevent  the 
foil  from  producing  excessive  pressure  induced  noise,  a  phenomenon  known  as  singing.  These 
modifications  were  made  based  on  recommendations  from  the  Naval  Surface  Warfare  Center, 
Carderock  Division  (NAVSEA),  who  developed  the  foil  design  for  the  University  of  Michigan. 
The  foil  model  used  for  the  project  had  a  chord  of  7  ft  and  a  span  of  10  ft.  The  foil  was  cast  from 
a  naval  grade  nickel-aluminum  bronze  and  machined  to  a  nominal  surface  roughness  of  25  pin. 
See  figure  1  for  a  diagram  of  the  foil  cross-section  and  figure  2  for  a  close  up  of  the  trailing  edge. 
Several  different  measurements  were  made  on  the  foil  over  the  course  of  the  experiment.  These 
included:  laser  doppler  velocimety  (LDV)  average  velocity  measurements  made  both  onboard 
the  foil  and  externally,  particle  imaging  velocimety  (PIV)  instantaneous  velocity  measurements 
of  the  trailing  edge  flow,  foil  vibration  levels  made  at  various  locations  in  the  foil  using  onboard 
mounted  accelerometers,  dynamic  pressure  fluctuation  measurements  made  at  the  trailing  edge 
of  the  foil  using  flush  mounted  dynamic  pressure  transducers  and  surface  static  pressure 
measurement  made  using  a  number  of  static  pressure  taps  and  differential  diaphragm  pressure 
transducers.  See  figure  3  for  an  overview  of  the  locations  of  the  various  onboard  instrumentation 
mounted  in  the  foil. 
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Figure  1:  Cross-Section  view  of  baseline  Hifoil 


Figure  2:  Detail  of  baseline  trailing  edge 


Leading  Edge 


Over  the  course  of  the  project  these  measurements  were  made  on  the  baseline  foil  at  flow 
speeds  ranging  from  0.25  m/s  through  18.3  m/s.  The  various  measurement  were  repeated  on  the 
baseline  foil  at  angles  of  attack  of  0  degrees,  -1  degrees  and  +1  degrees,  measured  relative  to  the 
flat  pressure  side  of  foil.  Additionally,  a  modification  was  made  to  the  trailing  edge  of  the 
baseline  foil  and  another  series  of  measurement  were  made  at  0  degrees  angle  of  attack  at  flow 
speeds  ranging  again  from  0.25  m/s  though  18.3  m/s.  See  figure  4  for  a  diagram  of  the  modified 
trailing  edge  relative  to  the  baseline  foil.  The  modification  is  represented  in  the  drawing  by  the 
alternate  hatch  pattern  seen  at  the  trailing  edge.  The  modification  consisted  of  a  piece  of  brass 
that  had  the  baseline  foil  shape  machined  into  its  underside  and  the  new  trailing  edge  shape 
machined  into  its  outer  surface.  The  modification  stretched  the  entire  span  of  the  foil  and  was 
attached  with  bolts  and  faired  into  the  original  trailing  edge  with  bonding  compound. 


Figure  4:  Detail  of  trailing  edge  modification 


The  Hifoil  project  was  conducted  in  three  phases.  Phase  one,  conducted  in  October  of 
1999  was  concerned  mainly  with  the  average  flow  field  at  the  trailing  edge  of  the  foil.  All 
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measurements  taken  during  this  phase  were  made  using  and  external  two-component  LDV 
system,  owned  and  operated  by  the  LCC  staff.  Phase  two  took  place  in  February  and  March  of 
2001  and  it  was  during  this  phase  that  the  majority  of  the  measurements  were  made  including 
PIV  at  the  trailing  edge,  foil  vibration,  onboard  LDV  measurements  of  the  boundary  layer  at  the 
mid-chord,  trailing  edge  pressure  fluctuations  and  surface  static  pressure  measurements.  It  was 
also  during  this  phase  that  the  angle  of  attack  study  was  performed.  Phase  three  was  conducted 
during  August  of  2001  and  was  concerned  with  making  low  speed  measurements  on  the  baseline 
foil  and  then  repeating  several  of  the  measurements  made  during  phase  two  on  the  foil  with  the 
modified  trailing  edge. 

Experimental  Facility 

All  three  phases  were  conducted  at  the  WBM-LCC  in  Memphis,  Tennessee.  The  LCC  is 
a  variable  pressure  recirculating  water  tunnel  with  can  flow  water  from  0.25  m/s  up  to  18.3  m/s 
through  at  43  ft  long  test  section  with  a  10  ft  by  10  ft  square  cross-section.  The  pressure  in  the 
test  section  can  be  varied  0.5  to  60  psia,  in  order  to  control  the  onset  of  cavitation  in  the  test 
section.  The  LCC  pump  is  powered  by  10,440  kW  motor.  Given  these  statistics,  a  foil  chord  of 
7ft  and  a  channel  water  temperature  range  of  70  to  104  °F,  the  chord  based  Reynolds  number 
possible  during  this  experiment  ranged  from  1  to  60  million,  using  the  definition: 

Rec=—  (1) 

v 

Where  U  is  the  free  stream  fluid  velocity,  c  is  the  chord  length  of  the  foil  and  u  is  the  kinematic 
viscosity  of  the  water,  which  is  a  function  of  water  temperature.  Figure  5  is  a  schematic  of  the 
LCC  and  figure  6  is  diagram  of  the  foil  model  mounted  in  the  test  section  of  the  LCC. 
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Figure  5:  Diagram  of  the  LCC 


FLOW 


CROSS  SECTION 


2)  All  dimensions  referenced  to  inside  shell  of  test  section 


Figure  6:  Drawing  of  the  Hifoil  mounted  in  the  LCC  test  section 


Pressure  Acquisition  System  Hardware 

In  order  to  measure  the  static  pressure  gradients  on  the  foil  it  was  necessary  to  design  and 
implement  and  entire  static  pressure  acquisition  system.  This  systems  consisted  of  the  following 
parts:  the  pressure  taps  themselves,  the  tubing  system  that  brings  the  fluid  out  of  the  test  section 
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of  the  channel  to  where  its  pressure  can  be  measured,  a  manifold  system  that  makes  it  possible  to 
valve  between  the  different  pressure  tap  lines  and  bleed  any  air  out  of  the  lines,  the  transducers 
used  to  make  the  pressure  measurements  and  finally  the  computer  system  which  records  the 
readings  from  the  transducers  and  controls  the  system. 

The  pressure  taps  are  the  most  important  aspect  of  this  system.  There  are  three  important 
considerations  to  be  made  in  deciding  how  to  measure  static  pressure  with  pressure  taps.  First, 
the  number  of  taps  to  be  used  is  critical.  The  surface  static  pressure  field  on  the  foil  is  obviously 
a  continuous  function  of  chordwise  position.  However,  it  is  impossible  to  measure  the  pressure 
continuously  because  this  would  take  an  infinite  number  of  pressure  taps.  Therefore  it  is 
important  to  decide  how  many  taps  will  adequately  capture  the  trends  of  the  pressure  field  while 
minimizing  the  number  of  taps  that  have  to  drilled  into  the  model  in  an  effort  to  minimize 
machining  costs.  In  preparation  for  the  Hifoil  project  a  simple  calculation  was  made  by 
NAVSEA  using  a  two-dimensional  representation  of  the  foil  in  the  LCC  test  section  to 
approximate  the  probable  pressure  field  on  the  surface  of  the  foil.  This  calculation  was  done 
using  simplified  boundary  conditions  and  a  two-dimensional  panel  method  integral  solver.  From 
the  calculation  it  was  determined  that  30  taps  would  be  required  to  accurately  capture  the  shape 
of  the  pressure  field.  A  larger  percentage  of  the  taps  had  to  be  located  near  to  the  nose  and  the 
trailing  edge  of  the  foil,  where  the  pressure  gradients  were  anticipated  to  be  the  largest.  Having 
more  taps  in  these  areas  will  make  it  easier  to  spatially  resolve  these  larger  gradients.  Figure  7 
shows  the  results  of  the  computation  made  by  NAVSEA  presented  as  the  negative  of  the 
nondimensional  pressure  coefficient,  Cp  versus  chordwise  position  in  percent  of  chord  (x/C). 

The  convention  is  to  plot  -CP  in  order  to  place  the  suction  side  data  curve  above  the  pressure  side 
curve  in  pressure  coefficient  plots  in  order  to  mirror  the  actual  geometry  of  the  foil.  Figure  7 
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also  show  the  x/C  location  of  each  of  the  30  taps,  while  table  1  shows  the  numeric  x/C  location 
of  each  tap.  The  dots  on  the  foil  cross-section  in  figure  7  indicate  the  location  of  the  pressure 
taps.  The  formula  used  for  computing  Cp  is  as  follows: 


C 


p 


P  -P 
- ~r  (2) 

kpVm~ 


Where  Pw  is  the  measured  wall  pressure,  Poo  is  the  free  stream  pressure,  p  is  the  fluid  density  and 
Uco  is  the  free  stream  velocity.  As  can  be  seen  from  this  equation,  the  pressure  coefficient  is  a 
nondimensional  way  to  represent  the  measured  wall  pressure  as  a  differential  from  the  free 


stream  pressure. 


Figure  7:  NAVSEA  predicted  pressure  distributions  and  pressure  tap  x/C  location 

It  was  also  important  to  carefully  determine  the  spanwise  location  of  each  tap  in  addition 
to  the  chordwise  location.  As  previously  stated,  putting  a  pressure  tap  in  a  model  to  measure  the 
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static  pressure  will  disturb  the  flow  in  that  spanwise  location.  Therefore,  no  piece  of 
instrumentation  can  lie  downstream  of  any  static  pressure  tap  without  the  upstream  pressure  tap 
producing  uncertainty  in  the  downstream  measurements.  This  made  it  necessary  to  place  each  of 
the  static  pressure  taps  in  a  slightly  different  spanwise  location  so  as  not  to  corrupt  any  other 
measurements.  Refer  back  to  figure  3  to  see  the  spanwise  location  of  each  pressure  tap  in  a 
planform  view. 

Lastly  the  size  of  each  pressure  tap  was  an  important  factor  to  consider.  The  correlation 
provided  by  Benedict  suggests  that  the  smaller  the  tap  size  the  better.  [3]  Smaller  pressure  taps 
disturb  the  surrounding  flow  less  and  result  in  smaller  hole  error  values.  However,  experience  in 
the  LCC  has  found  that  pressure  taps  less  than  one  sixteenth  of  an  inch  in  diameter  tend  to  clog 
with  the  random  debris  that  exists  in  the  water  used  in  the  LCC.  For  this  reason,  it  was  decided 
that  the  majority  of  the  pressure  taps  in  the  foil  should  be  one  sixteenth  of  an  inch  in  diameter  to 
avoid  clogging.  However,  near  the  leading  edge  of  the  foil  the  high  radius  of  curvature  of  the 
nose  required  a  slightly  smaller  size  tap  diameter  in  order  to  reduce  error.  A  one  sixteenth  inch 
diameter  tap  size  near  the  nose  of  the  foil  would  create  a  large  defect  in  the  surface  of  the  foil 
and  it  was  speculated  that  taps  of  this  size  would  most  likely  disturb  the  flow  too  much  at  the 
leading  edge  of  the  foil.  For  this  reason  is  was  determined  that  it  would  be  better  to  have  smaller 
taps  near  the  nose  and  risk  having  them  get  clogged  over  the  course  of  the  experiment  in  order  to 
reduce  the  possibility  of  corrupting  any  other  measurements  to  be  made  over  the  course  of  the 
project.  Table  1  also  shows  the  diameter  of  each  of  the  thirty  foil  static  pressure  taps. 
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Pressure  Side  Taps 


Tap 

x/C 

Diameter  (in) 

1 

0.0313 

2 

0.0313 

0.004 

0.0313 

4 

0.0313 

5 

0.0313 

6 

0.060 

0.0625 

mm 

0.0625 

8 

0.0625 

0.300 

0.0625 

10 

0.0625 

ii 

0.0625 

12 

0.800 

0.0625 

13 

0.880 

0.0625 

14 

0.930 

0.0625 

15 

0.958 

0.0625 

16 

0.975 

0.0625 

17 

0.988 

0.0625 

Tap 

X/C 

Diameter  (in) 

18 

0.985 

0.0625 

19 

0.975 

0.0625 

20 

0.958 

0.0625 

El 

m 

0.0625 

El 

mm 

0.0625 

0.500 

0.0625 

El 

0.300 

0.0625 

mm 

EKES 

0.0625 

■ 

0.0625 

El 

BBl!i 

0.0313 

IBflil 

0.0313 

E52I 

0.0313 

30 

0.001 

0.0313 

Table  1:  x/C  location  and  hole  size  of  foil  pressure  taps 


With  the  size  and  location  of  each  of  the  taps  determined,  the  next  challenge  was  to 
devise  a  method  for  actually  measuring  fluid  pressure  at  each  tap.  Within  the  foil,  the  opposite 
side  of  each  pressure  tap  was  machined  to  a  female  sixteenth  inch  national  pipe  taper  (NPT) 
thread.  A  special  adapter  fitting  was  used  to  mate  the  NPT  thread  to  an  o-ring  tube  fitting.  This 
made  it  possible  to  use  eighth  inch  outer  diameter  flexible  nylon  1 1  tubing  to  connect  to  each  of 
the  taps.  Eighth  inch  outer  diameter  tubing  was  chosen  to  minimize  the  amount  of  space 
occupied  by  static  pressure  tubing  within  the  foil.  Outside  of  the  test  section  of  the  water  tunnel 
each  of  the  thirty  tubes  coming  from  the  pressure  taps  was  then  connected  to  a  custom  built 
common  manifold. 

The  manifold  was  designed  with  two  main  purposes  in  mind.  First,  it  had  to  be  capable 
of  bleeding  the  air  out  of  each  of  the  pressure  tap  lines  and/or  capable  of  pushing  pressurized 
water  out  of  each  of  the  pressure  tap  lines  in  an  effort  to  clear  all  the  air  and  debris  from  each  of 


16 


the  pressure  lines.  This  was  necessary  because  the  presence  of  bubbles  or  other  debris  in  the 
pressure  tap  lines  could  have  corrupted  the  pressure  measurements.  Bubbles  in  the  fluid  lines 
could  cause  the  surface  pressures  to  be  inaccurately  transmitted  down  the  line  to  the  pressure 
transducers.  The  common  manifold  also  needed  to  be  able  to  isolate  each  pressure  tap  lines 
individually  in  order  to  make  it  possible  to  measure  the  pressure  in  each  individual  line 
independently.  There  are  basically  two  ways  to  do  accomplish  this  function.  Either  each  line 
can  have  its  own  dedicated  sensor  or  one  sensor  can  be  used,  but  it  must  then  be  possible  to 
isolate  each  pressure  line  from  the  sensor  with  valves.  The  first  method  was  not  cost  effective 
both  in  terms  of  number  of  transducers  that  would  be  required  for  the  thirty  foil  pressure  taps, 
nor  in  terms  of  the  time  it  would  take  to  calibrate  thirty  different  transducers  before  any  data 
could  be  collected.  For  this  reason  the  second  method  was  chosen.  As  a  result  the  common 
manifold  was  designed  and  built  with  the  capability  to  valve  on  or  off  each  of  the  pressure  lines 
independently.  The  process  of  isolating  each  fluid  line  with  valves  could  be  done  by  hand  with 
mechanical  Swagelok  compression  fitting  ball  valves.  Three  way  valves  were  used  to  direct  the 
fluid  pressure  in  one  of  two  directions,  either  to  a  common  bleed  manifold  which  could  be  sealed 
off  or  opened  to  atmospheric  pressure,  or  to  the  transducer  manifold. 

The  transducer  manifold  consisted  of  a  Scanivalve  solenoid  driven  rotary  valve.  This 
device  was  capable  of  taking  forty-eight  input  fluid  lines  and  opening  up  those  input  lines  to  a 
single  output  line  one  at  a  time.  This  device  could  be  controlled  with  digital  computer  logic, 
making  it  possible  to  computer  control  the  pressure  measurement  process.  Information  about  the 
Scanivalve  can  be  found  in  appendix  B. 

The  Scanivalve  system  was  used  to  make  it  possible  to  automatically  isolate  an  individual 
pressure  line,  a  necessary  capability  when  attempting  to  measure  static  pressure  on  several 
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pressure  taps  with  a  single  pressure  transducer.  The  Scanivalve  option  was  chosen  because  it 
was  more  cost  effective  than  building  two  separate  manifolds  and  an  individual  solenoid  driven 
ball  valve  for  each  pressure  tap  line,  a  system  that  would  have  cost  approximately  $400  per 
pressure  tap  and  required  a  entire  control  circuit  for  each  valve.  The  complete  Scanivalve  system 
cost  only  $2000  total  and  could  support  up  to  forty-eight  pressure  taps.  However,  the  Scanivalve 
device  proved  to  have  limited  reliability.  When  it  was  operational  the  Scanivalve  worked  very 
well,  but  it  did  tend  to  break  down,  get  clogged  with  debris  due  to  its  small  size  and  have 
electronic  control  problems  from  time  to  time. 

The  Scanivalve  system  uses  a  custom  0.063  inch  outer  diameter  tubing  in  all  of  its  fittings 
in  order  to  facilitate  the  geometric  constraints  of  the  system.  Custom  0.063  inch  to  quarter  inch 
outer  diameter  Swagelok  fittings  were  purchased  in  order  to  incorporate  the  Scanivalve  into  the 
pressure  acquisition  system.  Swagelok  quarter  inch  outer  diameter  tube  fittings  were  used  for  all 
the  other  the  plumbing  applications  in  the  static  pressure  manifold.  Swagelok  compression 
fittings  were  chosen  for  their  excellent  quality  and  reliability.  Brass  fittings  were  chosen  based 
on  cost.  This  choice  was  justified  because  brass  is  an  acceptable  material  for  use  in  low-pressure 
water  applications.  Quarter  inch  outer  diameter  tubing  was  chosen  due  to  the  greater  inner 
diameter  when  compared  to  the  eighth  inch  tubing  used  inside  the  foil.  The  greater  inner 
diameter  made  it  easier  to  push  out  trapped  bubbles  and  debris.  This  dimension  also  made  it  less 
likely  that  any  bubble  that  might  be  trapped  in  the  tubing  would  be  large  enough  to  entirely  block 
the  tube,  a  situation  that  could  create  errors  in  the  pressure  measurements.  A  special  reduction 
union  tube  fitting  was  used  to  step  up  each  of  the  eighth  inch  outer  diameter  tubes  coming  out  of 
the  foil  to  the  quarter  inch  outer  diameter  tubing  used  in  the  manifold. 


18 


Appendix  C  shows  a  schematic  of  the  pressure  acquisition  manifold  plumbing  system. 
The  system  could  support  up  to  forty-eight  pressure  taps  to  match  the  forty-eight  inputs  into  the 
Scanivalve.  It  was  designed  to  be  flexible  in  the  number  of  inputs  used  to  support  other  later 
experiments  that  would  perhaps  use  a  different  number  of  pressure  taps  than  the  Hifoil  project. 
The  manifold,  valves,  Scanivalve,  all  electronic  circuits  and  transducers  were  all  mounted  on  a 
self-contained  stand  built  from  Unistrut  framing  material,  sheet  metal  and  plywood.  Figure  8 
shows  a  photograph  of  the  pressure  acquisition  stand. 


Figure  8:  Photograph  of  static  pressure  acquisition  system  (D.  Bourgoyne) 
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For  the  second  phase  of  the  Hifoil  project  all  forty-eight  inputs  were  used,  including  the 
thirty  taps  on  the  model  and  an  additional  eighteen  wall  taps  to  record  the  wall  pressures  inside 
the  LCC  test  section  at  various  places  upstream  and  downstream  of  the  foil.  These  pressures 
were  recorded  to  provide  boundary  condition  information  for  the  Hifoil  project  data  set.  During 
phase  three  pressure  data  was  taken  on  only  the  thirty  foil  taps.  The  LCC  staff  collected  the  wall 
pressure  data  using  an  LCC  designed  and  operated  static  pressure  acquisition  system.  Wall 
pressure  data  will  not  be  addressed  in  this  document. 

To  actually  measure  the  pressures  in  each  tap,  a  differential  pressure  measurement 
approach  was  used.  In  the  case  of  the  Hifoil  project,  the  free  stream  pressure  was  taken  to  be  the 
pressure  measured  at  a  channel  wall  pressure  tap  that  lied  approximated  2  chord  lengths 
upstream  of  the  foil  in  the  LCC  test  section.  The  data  for  all  static  pressure  taps  was  taken  as  a 
differential  from  the  pressure  measured  at  this  upstream  wall  tap  using  one  of  several  differential 
pressure  instruments  employed  over  the  course  of  the  experiment.  In  other  words,  each 
transducer  had  two  inputs  that  it  measured  a  pressure  differential  between.  One  input  into  the 
transducers  came  from  the  output  of  the  Scanivalve  used  in  the  pressure  acquisition  manifold 
while  the  other  input  came  from  the  upstream  reference  tap. 

Before  selecting  any  transducers  for  the  Hifoil  project,  a  preliminary  uncertainty  analysis 
was  done  to  take  into  account  the  range  of  error  expected  in  the  measurements.  The 
uncertainties  included  in  the  calculations  included  the  transducer’s  manufacturer  stated 
uncertainty,  analog  to  digital  quantization  error  due  to  the  data  acquisition  computer  limitations 
and  the  expected  static  pressure  hole  error  in  the  measurement.  Each  separate  transducer 
considered  had  a  slightly  different  stated  uncertainty  associated  with  it  as  reported  by  the 
manufacturer.  This  uncertainty  is  reported  in  a  certain  percentage  of  the  full  scale  operating 
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range  of  the  transducer.  This  error  source  was  usually  the  largest  contribution  to  the  total 
predicted  error  for  any  given  measurement.  The  quantization  error  in  any  measurement  was  a 
function  of  the  analog  to  digital  converter  (ADC)  used  to  convert  the  analog  voltage  output  from 
the  transducer  into  the  digital  signal  that  is  stored  by  the  computer.  This  uncertainty  is  a  function 
of  the  resolution  of  ADC  and  the  full  scale  range  of  voltages  output  by  the  transducer.  The  data 
acquisition  computer  system  will  be  discussed  in  detail  later  in  this  document,  however  it  was 
found  that  the  quantization  error  in  any  measurement  was  orders  of  magnitude  less  than  the  other 
elemental  sources  of  error,  and  was  thus  essentially  not  a  factor  in  the  total  uncertainty  of  any 
measurement.  The  only  other  expected  source  of  error  in  the  static  pressure  measurements  was 
static  pressure  tap  hole  error.  Based  on  the  previous  experimental  data  discussed  earlier  and 
Benedict’s  findings,  hole  error  is  known  to  be  a  bias  error  that  is  a  function  of  hole  size  and 
shape,  flow  speed  and  fluid  properties.  [3]  Using  the  known  empirical  relations,  an  estimation  of 
static  pressure  hole  error  was  included  in  the  preliminary  uncertainty  analysis.  A  sample 
worksheet  of  this  preliminary  calculation  is  contained  in  table  2. 


U  (m/s) 

18.3 

12 

6 

3 

Transducer  Full  Scale 

-Rosemount  3051  CD 

Error  TvDes- 

hole  ( +/-  psi) 

0.1457 

0.0501 

0.0076 

0.0004 

transducer  (+/-  psi) 

0.0555 

0.0555 

0.0555 

0.0555 

-error  (mfg)  0.00075 
-FS  (psi)  74 

A/D  (+/-  psi) 

-Efsr(V)  8 

-QE  (+/-  V)  0.001953 

0.0156 

0.0156 

0.0156 

0.0156 

Total  (+/- psi) 

0.1566 

0.0763 

0.0581 

0.0576 

Table  2:  Sample  calculation  of  expected  uncertainty  in  pressure  measurements 
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Based  on  this  type  of  a  rough  calculation,  two  acceptable  candidate  transducers  were 
selected  due  to  their  specifications  and  ability  to  measure  differential  pressure  in  water.  The  two 
models  included  the  Setra  P230  differential  pressure  transducer  and  the  Rosemount  305 1CD 
differential  pressure  instrument.  Both  models  were  eventually  purchased  and  used  at  some  point 
over  the  course  of  the  project.  The  Setra  devices  were  originally  selected  due  to  their  more 
reasonable  cost.  However  the  Rosemount  instruments,  though  more  expensive,  proved  to  be 
much  better  measurement  devices.  Appendix  D  contains  pictures  and  general  specifications  for 
both  instruments. 

Selecting  a  single  pressure  transducer  for  adequately  measuring  the  foil  static  pressure 
over  the  entire  range  of  flow  speeds  of  the  test  proved  to  be  impractical.  The  largest  pressure 
differential  seen  by  the  foil  was  almost  25  psid,  which  is  the  stagnation  pressure  at  18.3  m/s, 
while  at  0.25  m/s  pressure  differentials  were  as  small  as  tenths  of  a  psi.  There  are  no 
commercially  available  differential  pressure  instruments  made  which  can  accurately  span  this 
entire  range.  This  made  it  necessary  to  use  a  number  of  transducers,  each  sized  for  a  certain 
range  of  pressures  and  intended  for  use  at  a  specific  set  of  flow  speeds.  For  phase  two  the 
anticipated  test  conditions  were  3,  6,  12  and  18.3  m/s,  thus  the  anticipated  pressure  range  was 
0.1-25  psid.  To  accommodate  this  range,  two  Setra  transducers  were  purchased,  one  with  a  +/- 
10  psid  range  for  the  smaller  pressure  differentials  and  a  +/-  50  psid  ranged  transducer  for  the 
higher  differentials.  A  solenoid  driven  ball  valve  was  used  to  isolate  the  lower  ranged  transducer 
from  the  line  pressure  when  the  line  pressure  was  higher  than  the  transducer’s  safe  operating 
range.  This  was  accomplished  by  closing  the  valve  after  each  measurement,  testing  the  line 
pressure  in  the  new  tap  line  with  the  larger  ranged  transducer  and  then  deciding  if  the  pressure 
differential  in  the  line  was  low  enough  to  not  damage  the  lower  ranged  transducer.  If  the 
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pressure  level  was  safe,  the  isolating  valve  was  opened  and  data  was  recorded  for  that  tap  on 
both  transducers,  while  if  the  pressure  was  too  high,  the  transducer  was  kept  isolated  and  data 
was  record  only  on  the  larger  transducer.  This  process  was  achieved  entirely  by  computer 
programming  and  control  and  is  discussed  later  in  the  document. 

A  single  Rosemount  transducer  was  also  purchased  for  phase  two  in  case  the  accuracy  of 
the  Setra  transducers  proved  to  be  unacceptable.  This  transducer  had  a  max  range  of  +/-  36  psid. 
The  Setra  transducers  cost  approximately  one  third  of  the  Rosemount  device  and  thus  were  a 
better  option  from  a  cost  perspective.  However,  the  Rosemount  transducers  were  known  to  be 
reliable  and  accurate  devices  and  so  one  was  purchased  as  a  backup.  As  could  have  been 
predicted,  the  Setra  devices  were  less  than  ideal  and  as  a  result  the  vast  majority  of  the  phase  two 
static  pressure  data  was  taken  on  the  single  Rosemount  device.  For  phase  three  an  additional 
Rosemount  instrument  with  a  max  operating  range  of  +/-  1  psid  was  purchased  to  more 
accurately  measure  the  static  pressure  differentials  at  the  lower  flow  speed  test  conditions  in  the 
phase  three  test  plan.  Efforts  were  also  made  to  devise  a  process  of  re-ranging  and  recalibrating 
the  Rosemount  transducers  to  reduce  their  full  scale  operating  ranges  and  the  resulting 
measurement  uncertainty.  The  Rosemount  device  has  the  capability  to  scale  down  its  operating 
range  using  a  special  computer  control  device  sold  by  Rosemount,  Inc.  This  allows  a  user  to  set 
the  range  of  the  transducer  to  whatever  is  needed  for  a  specific  measurement.  The  computer 
device  then  adjusts  the  transducer’s  internal  calibration  and  gain  level  to  accommodate  the 
desired  range.  For  instance,  if  the  expected  pressures  at  a  given  condition  predicted  pressures 
ranging  from  -5  up  to  12  psid  at  a  given  flow  speed,  the  +/-  36  psid  Rosemount  transducer  could 
have  its  range  spooled  down  to  +/- 13  psid,  making  the  full  scale  operating  range  26  psid,  which 
would  thus  reduce  the  uncertainty  in  the  measurements  made  at  that  test  condition.  This  process 
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was  used  throughout  phase  three,  and  table  3  shows  the  range  that  each  transducer  was  set  to  for 
each  flow  speed  in  phase  three  test  plan.  The  ranges  were  selected  to  minimize  the  error  in  the 
measurements  at  each,  test  condition.  This  approach  was  employed  in  an  effort  to  improve  the 
uncertainty  in  the  measurements  made  during  phase  three  over  those  made  during  phase  two 
when  the  transducer  ranges  were  the  same  throughout  the  test,  especially  at  the  low  speeds  where 
the  pressure  differentials  are  very  small. 


Transducer  Limits  (psi) 

U  (m/s)  Range  12 


<1 

A 

3.61 

0.09 

1.5 

B 

3.61 

0.36 

3 

C 

3.61 

0.90 

6 

C 

3.61 

0.90 

12 

E 

25.25 

0.90 

18.3 

E 

25.25 

0.90 

Table  3:  Phase  three  transducer  ranges  as  a  function  of  flow  speed 

Data  Acquisition  Computer  System  and  Software 

The  entire  pressure  acquisition  system  was  controlled  by  a  Windows  2000  based  personal 
computer  running  National  Instruments’  Labview  software.  A  National  Instruments  PCI-607  IE 
series  data  acquisition  board  was  used  to  interface  between  the  computer  and  the  hardware  in  the 
pressure  acquisition  system.  The  PCI-6071E  board  contains  a  12  bit  ADC,  32  differential  analog 
inputs  and  up  to  8  digital  input/output  (I/O)  channels  for  computer  control  of  various  processes. 
The  static  pressure  system  used  two  analog  inputs  to  record  data  from  the  two  transducers,  and  6 
digital  channels.  The  digital  channels  were  used  to  control  the  position  of  the  Scanivalve,  the 
isolation  valve  for  the  lower  range  pressure  transducer  and  to  read  the  position  of  the  Scanivalve 
using  an  encoder  built  into  the  Scanivalve  by  the  manufacturer. 

Control  and  data  acquisition  programs  were  written  for  the  static  pressure  system  in  the 
Labview  language.  Labview  is  a  National  Instruments  product  that  works  very  well  when  used 
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with  National  Instruments  hardware.  Labview  is  a  computer  programming  language  with  a 
graphical  interface  that  makes  is  very  easy  to  control  mechanical  processes,  acquire  data  from 
sensors,  perform  limited  data  processing,  view  data  in  near  real  time  and  save  data  to  files. 

There  were  three  major  programs  developed  for  the  Hifoil  project:  a  data  acquisition  program,  a 
transducer  calibration  program  and  a  data-sampling  program. 

The  first  program  controlled  the  acquisition  process  of  the  static  pressure  data.  This 
program  was  developed  to  take  a  single  set  of  data  for  a  given  test  condition  when  prompted  to 
do  so  by  the  user.  For  a  given  test  condition,  when  given  a  specific  number  of  pressure  taps  to 
take  data  from,  the  program  was  designed  to  automatically  select  the  appropriate  pressure  tap, 
take  a  time  series  of  data,  find  the  mean  and  standard  deviation  of  that  time  series,  record  that 
data  point  and  then  adjust  the  hardware  and  Scanivalve  for  the  next  tap  in  order  to  repeat  the 
process.  The  program  produces  a  data  file  that  contains  the  average  differential  pressure  and  the 
standard  deviation  of  the  time  series  as  measured  on  both  of  the  transducers,  for  each  desired 
pressure  tap  as  specified  at  the  beginning  of  the  run  by  the  user. 

Over  the  course  of  the  Hifoil  project,  three  separate  derivations  of  this  program  were 
developed.  In  the  original  program  the  computer  completely  controlled  the  data  acquisition 
process.  The  user  simply  input  the  flow  conditions,  the  sampling  rate,  number  of  samples 
desired  in  each  time  series  and  the  number  of  taps  on  which  to  take  data.  The  computer  than 
automatically  took  a  time  series  of  data  on  each  channel  sequentially,  found  the  mean  of  each 
time  series  and  saved  a  data  file  to  the  computer’s  hard  drive  that  contained  the  input  flow  data 
and  the  acquired  mean  differential  pressure  data.  The  first  derivation  of  this  program  was 
developed  to  allow  the  user  to  check  the  position  of  the  Scanivalve  before  allowing  the  computer 
to  take  data  on  that  channel.  Due  to  reliability  issues  with  the  Scanivalve  it  became  important  to 
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be  able  to  let  the  computer  take  data  and  advance  the  Scanivalve  to  the  next  channel,  but  then 
allow  the  user  to  verify  that  the  Scanivalve  position  had  actually  changed  to  the  correct  new 
position  before  allowing  the  computer  to  proceed  with  the  data  acquisition  process.  The  third 
derivation  of  the  program  was  written  to  allow  data  to  be  taken  in  the  event  that  the  Scanivalve 
was  inoperable.  In  this  event  the  Scanivalve  was  capped  off  and  the  transducer  manifold  was 
connected  to  the  common  bleed  manifold.  Then  the  mechanical  ball  valves  were  set  so  that  all  of 
the  pressure  tap  lines  except  for  one  were  fed  to  the  blocked  off  Scanivalve.  The  pressure  in  the 
manifold  was  thus  able  to  equalize  to  the  pressure  in  the  single  line  that  was  fed  into  the 
manifold  and  the  manifold  pressure  could  then  measured.  The  program  that  controlled  this  data 
acquisition  system  simply  took  and  recorded  data,  as  there  was  no  control  involved  in  this 
process.  The  program  would  take  the  data  on  a  channel  and  indicate  to  the  user  that  it  was  done. 
The  user  would  then  close  the  valve  of  the  completed  channel  and  open  the  valve  of  the  next 
channel  manually.  At  this  point  the  manifold  would  come  to  the  pressure  of  the  new  channel  and 
the  user  would  then  tell  the  computer  to  take  the  next  data  point.  The  net  result  of  this  program 
was  a  data  file  identical  to  the  data  files  produced  by  the  original  program,  however  this  program 
required  much  more  active  participation  by  the  user.  These  three  programs  were  written  to  allow 
for  a  variety  of  failures  that  could  occur  in  the  static  pressure  acquisition  system. 

The  second  program  developed  for  the  Hifoil  project  was  a  transducer  calibration 
program.  This  was  written  to  allow  the  user  to  calibrate  each  of  the  transducers  independently 
using  a  pressure  standard.  Both  types  of  transducers  stated  to  have  linear  calibrations  within 
their  full  scale  operating  range.  As  a  result,  these  calibrations  were  experimentally  verified 
based  on  the  linear  assumption  to  check  and  record  the  slope  and  zero  offset  for  all  transducers 
on  every  day  that  they  were  to  be  used  to  take  data.  This  linear  calibration  process  was  based  on 
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eleven  data  points,  5  with  the  known  pressure  applied  to  one  side  of  the  differential  pressure 
instrument,  five  with  a  known  pressure  applied  to  the  opposite  side  of  the  transducer  and  one 
zero  point  where  both  sides  of  the  transducer  were  simply  left  open  to  atmospheric  pressure,  thus 
resulting  in  a  zero  pressure  differential.  The  program  allowed  the  user  to  input  the  applied 
pressure  at  one  of  the  eleven  points.  The  computer  would  then  measure  and  record  the  average 
voltage  output  of  the  transducer  over  a  specified  period  of  time,  usually  3  seconds.  After  all 
eleven  points  are  recorded,  the  program  then  calculated  the  slope  and  intercept  of  the  calibration 
and  saved  the  raw  data  to  the  computer’s  hard  drive. 

The  last  program  written  was  a  simple  routine  that  plotted  the  real  time  pressure  in  a 
selected  line  as  a  function  of  time.  This  program  allowed  the  user  to  select  a  specific  pressure 
tap  and  then  see  the  pressure  in  that  line  as  measured  by  each  of  the  transducers.  This  program 
allowed  the  user  to  see  if  the  pressure  was  relatively  stable  in  time  and  to  see  how  the  calibration 
on  the  two  transducers  compared  for  a  given  pressure. 

The  programs  used  during  the  Hifoil  project  were  designed  to  make  data  acquisition 
easier  and  less  time  intensive  for  the  experimentalist,  however  the  ability  to  computer  control 
some  aspects  of  the  system  was  slightly  limited  by  hardware.  Specifically,  the  Scanivalve  was 
limited  in  how  it  could  move  between  its  forty-eight  inputs.  The  forty-eight  inputs  of  the 
Scanivalve  are  set  up  in  a  circular  pattern.  The  Scanivalve  system  is  set  up  to  index  the  valve 
from  one  input  to  the  next  input  sequentially  along  the  radius  of  the  circle  every  time  the 
Scanivalve  controller  sees  a  digital  logic  high  pulse  of  approximately  five  volts  as  produced  by  a 
computer  digital  output  signal.  One  of  the  inputs  was  randomly  selected  by  the  manufacturer  to 
be  the  number  one,  or  ‘home’  input,  and  the  rest  are  numbered  in  radial  order  from  the  number 
one  input.  The  result  of  this  is  that  it  is  a  relatively  easy  process  to  go  from  input  one  to  two  or 
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three,  etc.  However,  it  is  difficult  to  go  between  inputs  that  are  not  in  sequential  order  or  close  to 
each  other  in  number  order.  It  was  found  that  the  most  efficient  method  was  to  step  through  the 
inputs  sequentially  because  this  required  the  least  amount  of  pulses  to  be  issued  to  the  Scanivalve 
controller  by  the  computer.  As  a  result  the  program  was  set  up  to  record  data  from  the  pressure 
taps  in  a  sequential  manner.  The  user  could  select  to  take  data  on  less  than  the  total  number  of 
taps,  but  had  to  set  the  Scanivalve  up  to  start  on  a  specific  channel  manually.  For  the  purpose  of 
the  Hifoil  project,  the  thirty  pressure  taps  on  the  foil  were  set  up  as  inputs  one  through  thirty  on 
the  Scanivalve.  In  other  words,  when  running  automatically,  the  Scanivalve  stepped  through  the 
inputs  connected  to  foil  taps  first.  To  facilitate  this  methodology,  each  pressure  tap  was 
assigned  a  number  from  one  to  forty-eight,  and  it  was  this  number  that  was  used  by  the  computer 
program  to  associate  with  a  specific  pressure  tap.  For  example,  if  the  user  desired  to  look  at  the 
data  on  the  suction  side  of  the  foil  at  the  mid-chord,  he  would  have  to  set  the  Scanivalve 
manually  to  input  10,  the  number  of  the  mid-chord  tap,  and  then  start  the  data  acquisition 
program.  Because  of  these  constraints,  the  data  acquisition  system  was  limited  to  taking  data  on 
specified  number  of  taps  sequentially  with  the  only  variables  being  number  of  taps  and  starting 
tap.  The  graphical  code  and  user  interface  screens  for  each  of  the  programs  discussed  can  be 
found  in  appendix  E. 

Data  Regression  Techniques 

Calculating  and  plotting  the  nondimensional  pressure  coefficient  data  was  a  relatively 
simple  task.  The  voltages  output  by  the  transducer  were  converted  to  pressures  using  the 
corresponding  calibration  data  for  the  specific  transducers  and  ranges  used  to  collect  each 
specific  data  set  by  the  data  acquisition  program  as  the  data  set  was  being  recorded.  As  a  result 
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each  data  file  contained  mean  pressure  differential  information  as  a  function  of  both  instrument 
used  and  pressure  tap.  The  standard  deviation  of  the  time  series  record  at  each  tap  was  also 
saved  in  each  data  file  again  as  a  function  of  transducer  and  pressure  tap.  Lastly,  the  flow  speed 
and  water  temperature  for  each  test  condition  were  also  recorded  in  the  file  with  the  pressure 
data.  Using  this  data  and  equations  1  and  2,  the  Reynolds  number,  pressure  coefficient  and  the 
resulting  statistical  uncertainty  of  the  pressure  coefficient  value  for  each  data  point  could  be 
calculated.  Then,  the  Cp  data  for  each  tap  could  be  plotted  as  a  function  of  x/C  location  on  the 
foil.  Separate  Cp  curves  were  created  for  each  different  test  condition. 

Each  set  of  Cp  data  was  corrected  using  a  zero-flow  data  point  taken  each  test  day.  To 
take  this  data,  the  flow  in  LCC  was  brought  to  0  m/s  and  a  regular  static  pressure  data  run  was 
made.  Ideally,  the  differential  pressure  on  each  tap  at  zero-flow  should  be  zero,  however  this 
was  never  the  case.  The  zero-flow  condition  was  taken  at  some  point  during  each  testing  day  to 
record  the  inevitable  offsets  that  were  present  in  the  system  due  to  drifting  in  the  calibration 
offsets  of  the  transducers  or  from  any  slight  discrepancies  in  the  physical  system  due 
environmental  conditions  such  as  possible  debris  in  the  lines  or  temperature  changes. 

The  zero-flow  data  was  also  used  to  see  if  there  was  a  possible  clog  or  problem  with  any 
of  the  taps  or  inputs  into  the  Scanivalve.  If  the  zero-flow  data  suggested  a  problem  with  a  tap, 
that  problem  could  be  investigated  before  any  inaccurate  data  was  taken  on  that  tap. 

Uncertainty  in  the  pressure  measurements  was  calculated  in  terms  of  a  nondimensional 
offset  from  the  Cp  value.  This  uncertainty  has  both  a  precision  and  bias  component.  The 
statistical  component  comes  from  the  standard  deviation  of  the  data  time  series  as  calculated  by 
the  pressure  acquisition  program.  The  bias  component  is  the  results  of  the  static  pressure  tap 
hole  error  in  the  measurement.  To  compute  this  component  of  the  uncertainty  the  static  pressure 
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hole  error  associated  with  each  and  every  pressure  measurement  was  calculated  using  Benedict’s 
correlation.  [3]  The  standard  deviation  error  and  the  pressure  tap  hole  error  were  combined 
using  the  root  sum  squares  method  and  reported  as  the  total  uncertainty  associated  with  each 
corresponding  Cp  value. 

In  addition  to  calculating  the  pressure  coefficient  at  each  tap  measure  for  every  test 
condition,  the  pressure  lift  coefficient,  CL,  and  the  pressure  drag  coefficient,  CD,  for  each  test 
condition  were  also  calculated.  Each  set  of  data  recorded  was  nominally  comprised  of  the 
pressure  measured  at  each  of  the  thirty  taps  for  a  single  given  trailing  edge  configuration,  angle 
of  attack,  flow  speed  and  water  temperature.  Based  on  this  data,  each  Cp  curve  has  an  associated 
chordwise  Reynolds  number,  angle  of  attack  and  trailing  edge  configuration.  This  pressure 
curve  can  then  be  integrated  to  get  the  lift  and  drag  coefficients,  as  integrating  Cp  over  x/C  gives 
Cl  and  integrating  CP  over  y/C  yields  CD.  This  drag  coefficient  is  not  very  useful,  as  pressure 
drag  is  typically  a  very  small  component  of  the  total  drag  on  a  foil.  As  a  result  the  Cd  values  are 
not  reported  in  this  document.  The  calculations  used  to  determine  these  values  were  carried  out 
using  the  Matlab  programming  language,  and  can  be  found  in  appendix  F.  In  order  to  integrate 
the  Cp  curves,  a  simple  spline  function  in  Matlab  was  used  to  interpolate  the  CP  values  at  x/C  and 
y/C  values  between  the  known  data  points.  Then  using  a  Matlab  trapezoid  rule  numerical 
integrator,  the  CP  curves  were  integrated  to  yield  corresponding  lift  and  drag  values  for  each  data 
set  at  a  specific  test  condition. 

.The  uncertainty  in  these  integrated  quantities  was  determined  by  splining  and  integrating 
the  error  bars  associated  with  the  CP  curves.  An  overstated  integral  value  and  understated 
integral  value  were  calculated  for  each  CP  curve.  The  difference  between  these  two  integral 
quantities  was  then  found,  and  this  value  is  reported  as  the  uncertainty  in  the  corresponding  lift 
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quantity.  This  method  is  simple,  but  it  does  produce  a  basic  idea  of  the  uncertainty  associated 
with  the  calculated  lift  values  reported  in  this  document.  These  calculations  can  also  be  seen  in 
the  Matlab  code  enclosed  in  appendix  F. 
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Results 


The  static  pressure  data  collected  during  the  Hifoil  project  is  presented  here  in  form  of 
negative  CP  as  a  function  of  x/C  location  on  the  foil.  In  general,  light  colored  data  points 
indicate  low  flow  speeds,  while  darker  colored  data  points  represent  higher  speed  flows.  As  a 
convention,  the  four  configurations  of  the  foil  tested  during  the  Hifoil  project  are  differentiated 
between  in  these  plots  by  different  shaped  markers.  Square  markers  (■)  indicate  data  taken  on 
the  baseline  foil  at  0  degrees  angle  of  attack,  gradients  (▼)  indicate  -1  degrees  angle  of  attack 
data  on  the  baseline  foil,  deltas  (a)  represent  data  taken  at  +1  degrees  angle  of  attack  on  the 
baseline  foil  and  diamonds  (♦)  indicate  data  taken  at  0  degrees  angle  of  attack  with  the  modified 
trailing  edge. 

The  Hifoil  project  covered  flow  speeds  from  0.25  to  18.3  m/s.  However,  due  to 
limitations  in  the  transducers  and  hardware  used  in  the  static  pressure  system,  below  1.5  m/s  the 
pressure  head  in  the  system  was  too  small  to  produce  repeatable  static  pressure  data.  For  this 
reason  the  static  pressure  data  presented  here  covers  flow  speeds  from  1.5  to  18.3  m/s  and 
temperatures  from  75  to  103  °F,  resulting  in  a  Reynolds  number  range  of  3.9  to  57  million.  See 
the  following  table  for  exact  Rec  ranges  for  each  test  speed  for  which  static  pressure  data  is 
presented. 


Etc 

U  (m/s)  Average  Max  Min 


1.5 

4.2E+06 

4.3E+06 

3.9E+06 

3 

8.5E+06 

9.3E+06 

7.7E+06 

6 

1 .7E+07 

1 .8E+07 

1.5E+07 

12 

3.3E+07 

3.8E+07 

3.0E+07 

18.3 

5.1E+07 

5.7E+07 

4.8E+07 

Table  4:  Rec  ranges  for  static  pressure  data  at  each  test  speed 


The  first  plot  shown  here  contains  all  data  taken  at  0  degrees  AOA  on  the  foil  with  its 
original  trailing  edge.  As  can  be  seen  in  this  plot,  the  curves  for  most  of  the  test  conditions 
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collapse  within  experimental  uncertainty,  suggesting  little  to  no  Reynolds  number  dependence  in 
the  static  pressure  data.  However,  there  are  a  few  curves  that  are  obviously  offset  from  the  bulk 
of  the  data.  These  curves  are  the  data  taken  at  1.5  and  3  m/s,  and  were  collected  on  the  same  day 
of  the  test  using  two  separate  calibrations.  All  three  1.5  m/s  data  sets  overlay  with  each  other 
and  all  four  of  the  3  m/s  data  also  collapse  on  each  other.  Additionally  there  exists  static 
pressure  data  at  3  m/s  taken  during  phase  two  which  does  collapses  with  the  curves  from  the 
other  flow  speeds.  This  3  m/s  data  is  not  presented  here  due  to  its  relatively  high  associated 
uncertainty,  however  it  clearly  suggests  that  the  3  m/s  data  should  overlay  with  the  other  flow 
speeds.  This  suggests  that  a  constant  offset  error  existed  in  the  static  pressure  system  on  the  day 
that  these  1.5  and  3  m/s  data  sets  were  collected. 


x/C 

Figure  9:  Original  trailing  edge  configuration  at  0  degrees  AO  A,  all  data,  uncorrected 
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Due  to  the  constant  and  repeatable  nature  of  this  offset  error  and  the  fact  that  the  errors 
all  occurred  on  the  same  day,  this  offset  is  most  likely  the  result  of  an  error  made  in  calibration. 
As  a  result  this  error  was  corrected  using  the  following  method.  A  12  m/s  data  run  that  collapsed 
with  the  majority  of  the  other  data  sets  was  selected  as  an  arbitrary  baseline  data  set.  Then,  the 
offset  of  each  Cp  data  point  in  the  questionable  curves  from  the  baseline  Cp  curve  was  calculated 
at  each  corresponding  x/C  location.  The  individual  offsets  at  each  of  the  30  chordwise  locations 
were  then  averaged  together  for  each  of  the  seven  data  sets.  Each  of  these  resulting  single 
average  offsets  was  then  applied  as  a  constant  correction  to  each  of  the  corresponding 
questionable  data  sets.  The  next  figure  presents  the  results  of  this  correction  method  as  plotted 
against  the  remaining,  uncorrected  curves. 


x/C 

Figure  10:  Original  trailing  edge  configuration  at  0  degrees  AOA,  all  data,  corrected 
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As  can  be  seen  by  the  agreement  in  this  plot,  the  offsets  have  now  been  eliminated  and 
the  1.5  and  3  m/s  data  collapses  with  the  rest  of  the  data  sets  at  this  foil  configuration  within 
experimental  uncertainty. 

The  next  four  plots  presented  show  the  data  collected  for  the  original  trailing  edge  shape 
at  -1  and  +1  degree  AOA.  This  data  was  all  taken  during  the  angle  of  attack  study  done  during 
phase  two,  and  as  a  result  it  has  some  differences  from  the  other  data  presented  in  this  document. 
First,  the  lowest  test  speed  in  phase  two  was  3  m/s,  so  no  data  exists  for  -1  and  +1  degrees  AOA 
at  1.5  m/s.  Second,  the  transducers  used  in  phase  two  had  greater  uncertainties  associated  with 
them  and  as  a  result  the  uncertainty  in  the  CP  values  from  phase  two,  especially  at  the  lower  test 
speeds,  are  substantially  higher  than  those  seen  in  the  phase  three  data.  The  following  figures 
show  first,  all  the  data  collected  during  the  angle  of  attack  study  for  a  specific  geometry,  and 
then  the  same  data  plotted  again,  this  time  without  the  3  m/s  results.  This  is  done  to  show  two 
things.  First,  that  the  3  m/s  data,  despite  its  large  uncertainty,  does  collapse  with  the  data  taken 
at  higher  flow  speeds,  and  second,  that  when  the  data  taken  with  the  incorrectly  ranged 
transducers  is  not  considered,  the  remaining  data  has  a  very  low  associated  uncertainty.  This 
evidence  suggests  that  while  the  uncertainty  in  the  3  m/s  data  is  high,  the  true  error  in  these  data 
sets  is  more  than  likely  only  a  fraction  of  the  stated  maximum  uncertainty. 
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The  next  plot  presented  shows  the  data  collected  at  0  degrees  AOA  on  the  foil  with  the 
modified  trailing  edge.  The  data  curves  for  this  configuration  closely  resemble  the  0  degrees 
AOA  data  on  the  original  foil.  Further  comparisons  between  these  two  geometries  will  be  made 
later  in  this  document. 


x/C 

Figure  15:  Modified  trailing  edge  configuration  at  0  degrees  AOA,  all  data 

The  next  three  plots  show  how  the  actual  data  collected  compares  with  the  pressure 
curves  predicted  by  the  preliminary  NAVSEA  computations.  They  are  included  to  simply  prove 
that  the  static  pressure  measurement  made  during  the  Hifoil  project  produced  results  close  to 
what  was  predicted  by  the  simple  preliminary  calculations.  Further  computational  fluid 
dynamics  (CFD)  will  be  done  by  various  groups  in  the  near  future  to  more  accurately  model  and 
study  the  data  collected  by  the  Hifoil  project.  The  plots  presented  here  show  the  actual  data  for 
the  original  foil  shape  at  0,  -1  and  +1  degrees  angle  of  attack.  The  3  m/s  data  has  been  excluded 
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Figure  17:  Original  TE,  -1  degrees  AOA  data  vs  NAVSEA  prediction 
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Figure  18:  Original  TE,  +1  degrees  AOA  data  vs  NAVSEA  prediction 
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The  next  two  figures  present  data  from  the  baseline  foil  at  various  angles  of  attack.  The 
data  is  presented  at  both  high  and  low  flow  speeds  to  show  that  the  shape  of  the  curves  at  both 
extremes  of  the  velocity  range  studied  are  similar.  These  plots  show  the  obvious  result;  that  the 
pressure  differentials  between  the  suction  and  pressure  sides  of  the  foil  increase  with  angle  of 
attack. 


Figure  19:  Original  TE,  All  angles  of  attack  at  3  m/s 


Figure  20:  Original  TE,  All  angles  of  attack  at  18.3  m/s 


The  next  four  figures  show  the  effect  that  the  trailing  edge  modification  had  on  the  static 
pressure  distributions.  Data  from  the  two  trailing  edge  configurations  is  compared  at  1.5  and 
18.3  m/s  in  order  to  see  the  effect  of  different  flow  speeds  on  the  trailing  edge  comparison.  First 
all  the  data  is  presented,  and  then  a  close  up  view  of  the  trailing  edge  data  is  presented.  This  is 
done  because  the  most  appreciable  deviation  between  the  compared  data  sets  occurs  near  the 
trailing  edge,  as  might  be  expected.  However,  based  on  the  following  plots  no  obvious  trends  in 
the  data  can  be  seen  as  a  result  of  the  trailing  edge  modification.  The  only  change  seen  in  the 
data  is  a  greater  scatter  in  the  data  points  at  the  trailing  edge;  a  result  that  is  more  likely  due  to 
the  changes  made  to  the  geometry  of  the  pressure  taps  themselves  than  to  the  changes  in  the  flow 
field  caused  by  the  new  trailing  edge  geometry.  Based  on  the  data  in  these  plots,  no 
generalizations  can  be  made  about  the  effect  of  changing  the  trailing  edge  shape.  However,  it 
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should  be  pointed  out  that  the  modification  to  the  trailing  edge  of  the  foil  was  very  small  and 
perhaps  the  resulting  changes  to  the  static  pressure  field  were  too  small  to  be  accurately 
measured  by  the  pressure  transducers  used  in  this  project.  On  the  close  up  views  of  the  data  at 
the  trailing  edge  of  the  foil,  data  points  taken  on  the  original  foil  are  linked  with  a  dashed  line, 
while  modified  trailing  edge  data  is  represented  by  a  solid  line  to  make  any  differences  between 
the  two  configurations  easier  to  see. 
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Figure  21:  0  degrees  AOA,  original  vs  modified  trailing  edge  data  at  1.5  m/s 


43 


Figure  22:  Close  up  of  the  trailing  edge  data  in  the  previous  figure 
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Figure  23:  0  degrees  AOA,  original  vs  modified  trailing  edge  data  at  18.3  m/s 
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Figure  24:  Close  up  of  the  trailing  edge  data  in  the  previous  figure 

The  following  figure  shows  the  results  of  the  lift  calculations  made  by  integrating  the  CP 
curves.  The  four  lines  on  this  plot  are  linear  best  fits  of  the  CL  for  each  of  four  test  geometries  in 
the  Hifoil  project.  The  large  error  bars  seen  at  the  low  Reynolds  number  are  the  result  of  the 
large  uncertainty  in  the  3  m/s  CP  data  curves  at  -1  and  +1  degrees  AOA  discussed  earlier.  This 
plot  shows  that  there  is  a  slight  increase  in  lift  with  Reynolds  number  at  0  and  +1  degrees  AOA 
on  the  original  foil.  It  also  show  a  slight  decrease  in  Cl  with  Reynolds  number  for  the  -1  degree 
AOA  configuration.  However,  both  of  these  slight  trends  fall  entirely  within  the  error  bars  of  the 
calculated  lift  values  and  cannot  be  argued  with  any  degree  of  certainty.  The  Cl  for  the  foil  with 
the  modified  trailing  edge  shows  essentially  no  Reynolds  number  dependence. 
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Figure  25:  Calculated  lift  coefficient  vs  Reynolds  number  for  all  configurations 


The  calculated  uncertainty  levels  in  the  Hifoil  static  pressure  data  are  reported  below  in 
table  5  and  figure  26.  Table  5  shows  both  the  maximum  and  minimum  uncertainties  found  in  the 
complete  data  set,  both  in  the  nondimensional  Cp  form  and  also  as  a  differential  pressure  value  in 
psid.  Figure  26  shows  both  the  maximum  CP  and  maximum  pressure  uncertainty  levels  as  a 
function  of  Re<;.  This  figure  shows  that  the  maximum  CP  uncertainty  comes  from  the  low 
Reynolds  number  data,  most  likely  from  the  -1  and  +1  degrees  angle  of  attack  data  taken  during 
phase  two  with  the  inaccurately  ranged  transducers.  The  figure  also  shows  that  in  terms  of 
pressure  the  greatest  uncertainty  lies  at  the  highest  Reynolds  numbers,  but  that  this  high 
dimensional  uncertainty  converts  to  a  relative  small  nondimensional  value  when  normalized  by 
flow  speed.  However,  the  maximum  uncertainty  reported  in  table  5  may  be  misleading.  At  0 
degrees  AOA  all  static  pressure  curves  collapse  within  a  maximum  scatter  of  0.03  on  the  Cp  plot. 
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This  claim  is  also  supported  by  the  bulk  of  the  CP  uncertainty  data  reported  in  figure  26  lying 
below  a  value  of  0.04.  This  would  suggest  that  an  uncertainty  of  0.03  to  0.04  is  more  indicative 
of  the  true  error  in  the  static  pressure  measurements. 


Max 

Min 


+/■  Cp  +/•  psi 


0.1143 

0.4435 

0.0093 

0.0082 

Table  5:Maximum  and  minimum  uncertainties  in  the  Hifoil  static  pressure  data 
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Figure  26:  Maximum  uncertainty  values  as  a  function  of  Reynolds  number 
Additional  static  pressure  data  plots  from  the  Hifoil  project  can  be  found  in  appendix  A. 
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Conclusions 


The  static  pressure  data  collected  during  the  Hifoil  project  proved  to  be  of  good  quality. 
The  data  agreed  reasonably  well  with  the  preliminary  pressure  data  prediction  as  calculated  by 
NAVSEA,  which  proved  that  the  measurements  made  were  at  the  very  least  on  the  right  order  of 
magnitude.  The  collapse  of  the  data  taken  at  various  Rec  values  suggests  that  the  surface 
pressures  were  measured  well  and  that  the  measurement  and  calibration  techniques  were 
reasonably  repeatable.  Errant  data  was  corrected  when  possible  with  techniques  such  as  using 
zero-flow  data  or  data  taken  on  other  test  days  at  similar  flow  conditions  to  correct  offset  errors. 
Uncertainty  levels  were  high  for  some  data  sets,  but  on  the  whole  the  uncertainty  associated  with 
the  measurements  made  was  acceptable  when  the  optimal  pressure  instrument  for  a 
corresponding  flow  condition  was  used. 

No  significant  Reynolds  number  dependence  exists  in  the  static  pressure  data  collected 
over  the  Reynolds  number  range  of  the  Hifoil  project.  This  is  supported  by  the  collapse  of  the 
Cp  curves  taken  at  various  Reynolds  numbers  and  the  very  slight  trends  in  the  calculated  lift 
coefficients.  Additionally,  the  data  presented  shows  the  expected  increase  in  lift  with  increase  in 
angle  of  attack.  Lastly,  no  appreciable  dependence  in  the  pressure  data  was  observed  as  a  result 
of  the  change  in  trailing  edge  geometry. 

Based  on  the  data  collected,  the  static  pressure  system  designed  for  the  Hifoil  project 
functioned  well  and  produced  an  accurate  and  useful  high  Reynolds  number  static  pressure  data 
set,  fulfilling  one  of  the  major  goals  of  the  Hifoil  project  for  the  U.S.  Navy. 
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Appendix  A:  Additional  Data  Plots 

Contained  in  this  appendix  are  additional  annotated  plots  of  the  static  pressure  data  taken 
during  the  Hifoil  project.  Individual  plots  at  each  test  condition  for  which  static  pressure 
data  was  taken  are  presented  here. 
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Figure  1:  Original  trailing  edge,  O  degrees  AOA,  1.5  m/s 
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Figure  8:  Original  trailing  edge,  -1  degrees  AOA,  12  m/s 
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Figure  14:  Modified  trailing  edge,  0  degrees  AOA,  1.5  m/s 
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igure  15:  Modified  trailing  edge,  0  degrees  AOA,  3  m/s 
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Figure  16:  Modified  trailing  edge,  0  degrees  AOA,  6  m/s 
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Figure  17:  Modified  trailing  edge,  0  degrees  AOA,  12  m/s 


Appendix  B:  Scanivalve  Information 


Below  is  an  annotated  image  of  the  J-Type  Scanivalve  used  on  the  Hifoil  project.  This 
picture  comes  from  the  Scanivalve  Corporation  website  and  briefly  describes  how  the 
device  works. 

Px  multiple  pressure  Input 
lines  to  be  measured 

480M  pneumatic  connector, 
quick  Installation  of  Scanivalve 

Pressure  transducer 
Typical  port 
Elastic  duct 
Collector  hole 
Stator 
Rotor 

Thrust  Bearing 

Position  transmitter  (Encoder) 
electrically  Indicates  port 
being  measured 

Scanivalve  drive 
(Rot  shown) 

The  next  image  is  a  schematic  of  the  solenoid  driver  that  powers  the  J-Type  Scanivalve. 
Below  the  image  is  technical  information  about  the  drive  unit. 


This  stepper  drive  consists  of  a  12  step,  size  4  Ledex  solenoid  driving  a  4-1  gear  reducer. 
This  combination  makes  48  steps  per  revolution.  It  auto  homes  48  steps  in  two  seconds. 
Both  wafer  switches  rotate  at  solenoid  rpm. 


This  drive  requires  30  milliseconds  for  spring  return.  Use  Scanivalve  Model  CTLR2/S2- 
S6,  which  is  a  solid  state  pulse  length  feed  back  driving  circuit  for  providing  24Vdc  to 
operate  the  solenoid  drive. 

ScanCo#  JS4-48 - 48  steps/revolution,  wired  for  auto  home 

SPECIFICATIONS 

Duty  Cycle:  1/4  at  24  Vdc  or  44  seconds  "on"  period. 

Resistance:  Available  with  8.03  ohm  coil  for  operation  at  24  Vdc. 

Switch  Decks:  One  auto  home  wafer.  Wafer  switches  have  solid  silver  contacts 
(lubricate  with  light  grease.) 

Temperature:  -60  degrees  C  to  125  degrees  C,  can  be  operated  up  to  175  degrees  C 
(at  solenoid  not  ambient) 

(From  www.scaniva1ve.com) 


Static  Pressure  System  Schematic 


Appendix  C:  Experimental  Setup  Schematics 


Appendix  D:  Transducer  Specifications 

Setra  Transducer 


Model  230  Specifications  || 

Electrical 

Output 

Voltage 

0  to  5  VDC 

0  to  10  VDC 

Electrical 

Output 

Current 

4  to  20  mA 

Accuracy 

(RSS 

Method) 

±0.25%  Full  Scale 

Type  of 
Pressure 

Differential  (can  be  wet  both  sides) 

Pressure 

Ranges 

Unidirectional: 

Oto  1,2,5, 10,  25,50,  lOOpsid 

Bi-directional: 

0  to  ±0.5,  ±1,  ±2.5,  ±5,  ±10,  ±25,  ±50  psid 

Thermal 

Effects 

Thermal  Compensation  °F  (°C): 

30  to  150  (-1  to  65) 

%FS/100°F(100°C)max.zero: 

±2.0  (±3.6) 

%FS/1 00°F(1 00°C)max.span: 

±2.0  (±3. 6) 

Media 

Gases  or  liquids  compatible  with  300  series  and  17-4PH 
stainless  steel,  and  Viton®  "O"  ring.  (Hydrogen  not 
recommended  for  use  with  17-4  PH  stainless  steel.) 

*Viton®  is  a  registered  trademark  of  DuPont  Dow  Elastomers 

(From  www.setra.comt 


Rosemount  Transducer 


3051-  MODEL  305 1C  DIFFERENTIAL  PRESSURE  TRANSMITTER 

•  Superior  performance:  0.075%  accuracy,  100:1  rangeability 

•  Differential  pressure:  calibrated  spans  from  0.1  inH  2  O  to  2000  psi  (0,25  mbar 
to  138  bar) 

•  SST,  Hastelloy  C  ®  ,  Monel  ®  ,  Tantalum  (CD,  CG  only),  or  gold-plated  Monel 
process  isolators 

•  Compact,  rugged,  and  lightweight  design  for  easy  installation 

•  Compound  Range  allows  negative  calibrations 


(From  www.rosemount.com') 


Appendix  E:  Labview  Code 


Attached  are  Labview  printouts  of  the  Labview  programs  developed  for  use  with  the 
static  pressure  system  on  the  Hifoil  project.  Both  the  user  interface  and  the  block 
diagram  for  each  program  have  been  included.  The  programs  have  been  attached  in  the 
following  order: 


1.  Static  Pressure  Menu.vi 

2.  Pressure  Taps  (Mod  6a).vi 

3.  Pressure  Taps  (Mod  6b).vi 

4.  Pressure  Taps  (Mod  6c). vi 

5.  Calibration  (mod  3).vi 

6.  Static  Pressure  Sampling.vi 


Menu  progam 

Scanivalve  controlled  data  acquisition 
User  monitored,  Scanivalve  controlled 
User  controlled  data  acquisition 
Transducer  calibration  program 
Real  time  data  sampling  program 


Static  Pressure  Menu.vi 
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Last  modified  on  11/5/2001  at  11:51  AM 
Printed  on  11/5/2001  at  11:53  AM 


Static  Pressure  Menu.vi 

C:\HIFOIL\Static  Press  Taps\Phase  III\VIs\Static  Pressure  Menu.vi 
Last  modified  on  11/5/2001  at  11:51  AM 
Printed  on  11/5/2001  at  11:53  AM 


Static  Pressure  Menu.vi 

C:\HIFOIL\Static  Press  Taps\Phase  III\VIs\Static  Pressure  Menu.vi 
Last  modified  on  11/5/2001  at  11:51  AM 
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pressure  Taps  (Mod  6a).vi 
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Appendix  F:  Matlab  Code 


Attached  is  a  printout  of  the  Matlab  m-file  entitled  SPprocessing.m.  This  program  was 
developed  to  process,  correct  and  plot  the  static  pressure  data  collected  during  the  Hifoil 
project.  All  data  contained  in  this  document  was  processed  using  this  code.  This 
program  used  a  log  file  entitled  SP_log.txt  to  keep  track  of  the  entire  set  of  static  pressure 
data  files.  This  log  file  is  located  in  appendix  G 


%add  needed  file  paths 

addpath( 'C:\JHamel\Hifoil\Phase  III\SP  Post  Processing  Fall  2001'); 
addpath( 'C:\JHamel\Hifoil\Phase  III\SP  Post  Processing  Fall  2001\Data'); 
addpath( 'C:\JHamel\Hifoil\Phase  III\SP  Post  Processing  Fall  2001\Data\No 
F low  Data ' ) ; 

%addpath(' C:\HIFOIL\Static  Press  Taps\SP  Post  Processing  Fall  2001'); 
%addpath( 'C:\HIFOIL\Static  Press  Taps\SP  Post  Processing  Fall  2001 
\Data ' ) ; 

%addpath( 'C:\HIFOIL\Static  Press  Taps\SP  Post  Processing  Fall  2001\Data 
\No  Flow  Data ' ) ; 


%file  path  for  resulting  data  files 

rpath=' C:\JHamel\Hifoil\Phase  III\SP  Post  Processing  Fall  2001\Results 

\  1  ; 

%rpath= ’ C : \HIFOIL\Static  Press  Taps\SP  Post  Processing  Fall  2001\Results 
\\- 


%specify  number  format 
format  long  g 

%set  dx/c  step  size  for  all  calcs 
s  =  0 . 0  0  5  ;  %  0 . 0  0  0 1  ; 

%read  log  file 

[inf ile,U, aoa, te, errl , err2 , rho, nu, Re , noflow, outf ilel, outf ile2 ,  runid,  ff ,  o 
Id, cal , T) =textread (' SP_Log.txt %s  %f  %d  %d  %f  %f  %f  %f  %f  %s  %s  %s  %s 
%f  %d  %d  %f ',' headerlines ', 1) ; 


%read  x/c  and  tap  data 

[xcO , ycO , xcnl , ycnl , xcpl , ycpl , hd, bcpO , bcp02 , bcpnl , bcppl] =textread 
px_c.txt',  '%f  %f  %f  %f  %f  %f  %f  %f  %f  %f  %f 'headerlines ' , 2) ; 

%create  output  file  for  cld  data 

cldf ilel=strcat (rpath, ' cldo . txt ' ) ; 

f idl=f open (cldf ilel ,  '  w' )  ; 

fprintf (f idl , ' Variables="Re_C"\r\n ' ) ; 

fprintf (f idl , ' "TE»\r\n' ) ; 

fprintf (f idl , ' "AOA"\r\n' ) ; 

fprintf (f idl, ' "C_L"\r\n' ) ; 

fprintf (f idl , ' "err~C_L"\r\n ' ) ; 

fprintf (f idl , ' "C_D"\r\n' ) ; 

fprintf (f idl , ' "err-C_D"\r\n’ ) ; 

fprintf (fidl, 1 Zone\r\n' ) ; 

fprintf (fidl, 'T="TE  0  AOA  0Ao"\r\n'); 

cldf ile2=strcat (rpath, ’ cld0m.txt ’ ) ; 

f id2=fopen (cldf ile2 , ’ w’ ) ; 

fprintf (fid2 , ' Variables= "Re_C"\r\n ' ) ; 

fprintf (fid2 , ' "TE"\r\n’ ) ; 

fprintf (fid2, '"AOA"\r\n ' ) ; 

fprintf (f id2 , ' "C_L"\r\n' ) ; 

fprintf (f id2 , ' "err-C_L"\r\n' ) ; 

fprintf (fid2, ' "C_D"\r\n' ) ; 

fprintf (f id2 , ’ "err-C_D"\r\n' ) ; 

fprintf (f id2 , ' Zone\r\n ' )  ; 

fprintf  (fid2, 'T=''TE  2  AOA  0*o"\r\n ' )  ; 

cldf ile3=strcat (rpath, ' cldnl . txt ' ) ; 
f id3=fopen (cldf ile3 , ' w' ) ; 
fprintf (fid3, 'Variables="Re_C"\r\n' ) ; 
fprintf (fid3,  "'TE"\r\n' ) ; 


fprintf (f id3 , 
fprintf (f id3 , 
fprintf (f id3 , 
fprintf (f id3 , 
fprintf (f id3 , 
fprintf (f id3 , 
fprintf (f id3 , 


’  " AOA" \r\n ' )  ; 

’  " C__L"  \  r\n  '  )  ; 

'  "err-C_L,f\r\n  1  )  ; 

'  "C_JD"\r\n '  )  ; 

1  "err-C_D "\r\n ’ ) ; 

'  Zone\r\n 1 ) ; 

1  T- "TE  0  AOA  -lAo"\r\n' )  ; 


cldf ile4=strcat (rpath, ' cldpl . txt ’ ) ; 

f  id4  =  f  open  (cldf  ile4  ,'w’); 

fprintf (fid4 ,  1  Variables^ nRe_Cn \r\n  ‘ )  ; 

fprintf  (fid4,  1  "TE»\r\nf )  ; 

fprintf (fid4, ' "AOA"\r\n ' ) ; 

fprintf (fid4#  ' "C_L»\r\n 1 ) ; 

fprintf  (f  id4  ,  1  ”err-C_L,,\r\n  1  )  ; 

fprintf ( f id4 , • "C_D"\r\n 1 ) ; 

fprintf  (fid4,  '  ’'err-C_D,'\r\n,  )  ; 

fprintf (f id4 ,  1 Zone\r\n '  )  ; 

fprintf  (fid4,  ‘^"TE  0  AOA  +lAo" \r\n ’  )  ; 


deltaf ile=s treat (rpath,  1  delta .pit ' ) ; 
f id5=fopen (deltaf ile, 1 w* ) ; 
fprintf (fids, ' Variables^ M i M\r\n ' ) ; 
fprintf (fidS, ' MUn\r\nI ) ; 
fprintf (fidS, ' MAOAM\r\n ’ ) ; 
fprintf  (fid5,  '  ,,TEn\r\nl  )  ; 
fprintf (fidS, ' ”Tn\r\n' ) ; 
fprintf (fid5 , 1 HRe_Cu\r\n ’ ) ; 
fprintf  (fidS,  '  ,,Calibrationl,\r\n  1  )  ; 
fprintf (fidS, 1 "dC_P"\r\n ‘ ) ; 
fprintf  (fidS,  ’  ,,dP,,\r\n*  )  ; 


errf ile=s treat (rpath,  1  err. txt 1 ) ; 
fid6=f open (errf ile, 'w1 ) ; 
fprintf (fid6, ' Variables^ n i " \r\n 1 ) ; 
fprintf  (fid6,  '  ,,Ul,\r\n'  )  ; 
fprintf (fid6, ’ "AOAM\r\n 1 ) ; 
fprintf  ( f  id6  ,  1  ,fTE,T\r\nJ  )  ; 
fprintf (fid6, ' "Re_C" \r\n * ) ; 
fprintf (f id6 , ' "Max  C_Pn\r\n* ) ; 
fprintf (fide, ' "Max  Pn\r\n'); 
fprintf (fid6, ' "Max  C_P-err"\r\n' ) ; 
fprintf (fide , ' "Max  P-err"\r\n ' ) ; 

clear  cldfilel  cldfile2  cldfile3  cldfile4  deltafile  errfile 


%main  computation  loop 
[dl, d2] =size (infile) ; 
for  i=l:dl 

clear  infilei  noflowi  outfileli  outfile2i  runidi  tap  tdl  stdl  td2 
std2  valve  ntap  ntdl  nstdl  ntd2  nstd2 

clear  nvalve  offsetl  offset2  cpl  cp2  tau  Rd  he  cperrl  cperr2  cpstdl 
cpstd2  ebl  eb2  discrete  cp  eb 
%load  data  file 
inf ile i= char (infile (i) ) ; 
noflowi =char (nof low (i) ) ; 

outf ileli-strcat (rpath, char (outf ilel (i) ) ) ; 
outf ile2i  =  strcat (rpath, char (outf ile2 (i) )  )  ; 
runidi  =  char (runid (i) )  ; 


if  old (i) =-l 


[tap, tdl , stdl , td2 , std2 ,  valve] = text read (inf i lei ,  1  %f  %f  %f  %f  %f  % 
d ' , ' header lines ' , 10) ; 
else 

[tap, tdl , stdl , td2 , std2 , valve] = text read (inf i lei , 1 %f  %f  %f  %f  %f  % 
f 1 ,  ' headerlines  1 , 11)  ; 
end 

%correct  pressures  with  noflow  file 
if  old ( i ) ==1 

[ntap , ntdl , ns tdl , ntd2 , nstd2 , nvalve] =t ext read (nof lowi , ' %f  %f  %f  %f 
%f  %d 1 ,  1  headerlines  1 , 10) ; 
else 

[ntap , ntdl, nstdl /ntd2,nstd2, nvalve] =textread (nof lowi,  *%f  %f  %f  %f 
%f  %f ' ,  1  header lines  1 , 11)  ; 
end 

of f setl=mean (ntdl) ; 
of f set2=mean (ntd2) ; 
tdl=tdl-of f setl ; 
td2=td2-of f set2 ; 

%convert  pressures  to  Cp 

cpl= (tdl* (101327/14.7))/ (0 .5*rho (i) *U(i) A2) ; 

Cp2= ( td2* (101327/14.7))/ (0 . 5*rho (i) *U(i) A2) ; 

%hole  error  calculations 
tau=0 . 5*rho (i)*U(i)^2*0.002; 

Rd= (sqrt (tau/ rho (i) ) *hd) /nu (i)  ; 
if  Rd<  =  3 85 

he= (tau*0 . 000157*Rd . "l . 604 ) / (0 . 5*rho (i) *U(i) A2) ; 
else 

he= (tau*0.269*Rd. A 0.353) / (0 .5*rho(i) *U(i) A2) ; 

end 

%compute  error  bars 

cperrl= (errl (i) * (101327/14.7) ) / (0 .5*rho (i) *U(i) A2) / 
cperr2= (err2 (i)* (101327/14.7))/ (0. 5* rho (i) *U (i) A2 ) ; 
cpstdl= (stdl* (101327/14.7) )/ (0.5*rho(i) *U(i) A2) ; 
cpstd2= (std2  * (1013  27/14.7))/ (0 .5*rho (i) *U(i) A2) ; 
ebl=sqrt (cpstdl . A2+cperrlA2+he . A2) ; 
eb2=sqrt (cpstd2 . x2+cperr2x2+he . *2) ; 
for  j  =1 : 30 

if  valve ( j ) ==0 
cp2  ( j ) =0; 
eb2  (j ) =0; 
cp  ( j  )  =cpl  (j  )  ; 
eb  ( j  )  =ebl  (j  )  ; 
else 

cp ( j ) =cp2  ( j )  ; 
eb ( j ) =eb2  ( j )  ; 

end 

end 

%fudge  factor  correction 
%cp=cp+ff (i) ; 

if  aoa ( i) ==0 
if  te (i) ==0 

dcp=-cp-bcp0 ' ; 
else 

dcp=-cp-bcp02  '  ; 

end 

elseif  aoa(i)==-l 
dcp=-cp-bcpnl ' ; 
else 

dcp=-cp~bcppl ' ; 


end 

cpcorr=mean (dcp) ; 
cpc=cp+cpcorr  ; 

dpcorr=cpcorr* 0 . 5*rho (i) *U (i) A2*14 .7/101327; 

dout  = [ i ; U ( i ) ; aoa ( i ) ;te(i) ;T(i) ;Re(i) ;cal(i) ;cpcorr' ; dpcorr ' ] ; 
fprintf (fid5, ' %i  %4.2f  %d  %d  %5.2f  %8.1f  %d  %6.5f  %6 . 5f \r\n ' , dout) 

%output  discrete  Cp  points 

discrete=  [cpl 1 ;ebl ' ;cp2 ' ; eb2 1 ; valve  1 ]  * ; 

%write  Cp  data  file 

cpdata= [xcO ' ;yc0 1 ; -cp ; -cpc ; eb] ; 

fid=fopen (outfileli, 1  w1  ) ; 

fprintf  (fid,  1 Variables=ux/c”\r\n '  )  ; 

fprintf  (fid,  1  ,,y/c"\r\n  '  )  ; 

fprintf (fid,  1 "-Cp"\r\n1  )  ; 

fprintf (fid, 1 n-Cp  corr"\r\n' ) ; 

fprintf (fid, 1 nErr"\r\n ' ) ; 

fprintf (fid, 1 Zone\r\n 1 ) ; 

fprintf (fid, 'T=M 1 ) ; 

fprintf (fid, runidi) ; 

fprintf (fid,  ' n\r\n'  )  ; 

fprintf (fid, '%5.4f  %5.4f  %6.5f  %6.5f  %6.5f  \r\n ’ , cpdata) ; 
f close (fid) ; 

maxerr=max (eb) ; 

maxerrp=maxerr*0 . 5*rho (i) *U(i) ^2*14.7/101327; 
maxcp=max (cp) ; 

maxp=maxcp*0 . 5*rho (i) *U (i) ^2*14 . 7/101327 ; 

errout= [i;U (i) ; aoa (i) ; te (i) ;Re (i)  ;maxcp 1 ;maxp  * ;maxerr 1 ;maxerrp 1 ] ; 
fprintf (fid6, ' %i  %4.2f  %d  %d  %8.1f  %6.4f  %6.4f  %6.4f  %6.4f\r 
\n 1 , errout ) ; 

%decide  AoA  of  data  file,  spline  Cp  data,  compute  CL 
if  aoa  (i) =  =  0 

%0  degrees  AoA 
xt=xc0 (1) : s :xc0 (17) ; 
xb=xc0 (30) : s :xc0 (18) ; 
yc  =  0 : S : 0 . 0 8 0 02  ; 

XCt=xc0 ([1:17]); 
xcb=xc0 (  [18 : 30] )  ; 
cpt  =  Cp  (  [1 : 17]  )  ; 
errt=eb ([1:17]); 
cpt_high=cpt - errt ; 
cpt_low=cpt+errt ; 
cpb=cp ([18:30]); 
errb=eb (  [18 : 30] )  ; 
cpb_high=cpb+errb ; 
cpb_l ow= cpb - e  r rb ; 

%spline  data 

f cpt=interpl (xct , cpt , xt ,  ' spline  1 )  ; 
fcpt_high=interpl (xct , cpt_high,xt , ’spline1) ; 
f cpt_low=interpl (xct , cpt_low, xt ,  1  spline  1 ) ; 
f cpb  =  interpl (xcb, cpb, xb,  ' linear '  )  ; 
f cpb_high=interpl (xcb, cpb_high,xb, ’linear1) ; 
f cpb_low=interpl (xcb, cpb^low, xb, ’linear1) ; 

%calculate  CL  and  error  bars 

cl=s* (-trapz (fcpt) +trapz (fcpb) ) ; 

cl_high=s* (-trapz (fcpt_high) +trapz (fcpb_high) ) ; 


cl_low-s* (-trapz (fcpt_low) +trapz (fcpb_low) ) ; 

elh=cl_high-cl; 

ell=cl-cl_low; 

ecl= (elh+ell) /2 ; 

%drag  coefficient  calculations 
%create  variables 
ycf =ycO  (  [1:10,25:30]  )  ; 
ycr=yc0 ( [11 : 17]  )  ; 
cpf  =  cp {  [1:10,25:30]  )  / 
errf =eb ( [1:10,25:30]  )  ; 
cpf_high=cpf +errf  ; 
cpf_low=cpf -errf  ; 
cpr=cp ([11:17])/ 
errr=eb ([11:17]); 
cpr_high=cpr-errr ; 
cpr_low=cpr+errr ; 

%spline  data 

f cpf =interpl (ycf , cpf , yc ,  ' spline  1 )  ; 

f cpf_high=interpl (ycf , cpf_high, yc ,  1  spline  1 ) ; 

f cpf_low=interpl (ycf , cpf_low, yc ,  *  spline ' ) ; 

fcpr=interpl (ycr, cpr,yc,  ' spline  1 ) ; 

f cpr_high=interpl (ycr , cpr_high, yc , ’spline1) ; 

fcpr_low=interpl (ycr, cpr_low,yc, 'spline') ; 

%calculated  Cd  and  error  bars 

cd=s* (trapz (cpf) -trapz (cpr) ) ; 

cd_high=s* (trapz (cpf_high) -trapz (cpr_high) ) ; 

cd_low=s* (trapz (cpf_low) -trapz (cpr_low) ) ; 

edh=cd_high-cd; 

edl=cd-cd_low; 

ecd= (edh+edl) /2; 

cldoutput= [te (i)  ‘  , aoa (i) 1 , cl ' , eel ' , cd ' , ecd ' ] ; 
elseif  aoa  (i) ==-l 
%-l  degrees  AoA 
xt=xcnl (2) : s :xcnl (17) ; 
xb=xcnl (1) : s :xcnl (18) ; 
yc=0.00274:S:0. 09900; 

%lift  coefficient  calculations 

%create  variables 

xct=xcnl ([2:17]); 

xcb=xcnl ( [18:30,1]  )  ; 

cpt =cp ([2:17]); 

errt=eb ([2:17]); 

cpt_high=cpt-errt ; 

cpt_low=cpt+errt ; 

cpb=cp ([18:30,1]); 

errb=eb (  [18 : 3  0 , 1  ]  )  ; 

cpb__high=cpb+errb; 

cpb_low=cpb-errb; 

%spline  data 

f cpt=interpl (xct , cpt , xt ,  ' linear 1 )  ; 
f cpt_high=interpl (xct , cpt_high, xt , 'linear') ; 
f cpt_low=interpl (xct, cpt_low,xt ,  1  linear ’ ) ; 
f cpb=interpl (xeb, cpb,xb, 'linear') ; 
f cpb_high=interpl (xeb, cpb_high,xb, 'linear1) ; 
f cpb_low=interpl (xeb, cpb_low,xb, 'linear') ; 
%calculate  CL  and  error  bars 
cl=s* ( -trapz (f cpt) +trapz (f cpb) ) ; 
cl_high=s*  (  -trapz  (f cpt_high)  H-trapz  (f  cpb_high)  ) 
cl_low=s* ( -trapz (f cpt_low) +trapz (f cpb_low) ) ; 
elh=cl_high-cl ; 
ell=cl -cl_low; 


ecl= (elh+ell ) /2 ; 

%drag  coefficient  calculations 

%create  variables 

ycf =ycnl ( [1:10,27:30]  )  ; 

ycr=ycnl ( [11:26]  )  ; 

cpf  =  cp (  [1 : 10,27:30]  )  ; 

errf =eb ( [1:10,27:30]  )  ; 

cpf_high=cpf +errf ; 

cpf_low=cpf -errf ; 

cpr=cp ( [11 :  26]  )  ; 

errr=eb {[11:26]); 

cpr_high=cpr-errr ; 

cpr_low=cpr+errr; 

%spline  data 

f cpf =interpl (ycf , cpf , yc , * spline ' ) ; 

f cpf_high=interpl (ycf , cpf_high, yc , ’spline') ; 

f cpf_low^interpl (ycf, cpf_low,yc, ’ spline ’ ) ; 

fcpr=interpl (ycr, cpr,yc, ’ spline ’ ) ; 

f cpr_high=interpl (ycr , cpr_high, yc , ’spline*) ; 

fcpr_low=interpl (ycr, cpr_low,yc, ’spline’) ; 

%calculated  Cd  and  error  bars 

cd=s* (trapz (cpf) -trapz (cpr) ) ; 

cd_high=s* (trapz (cpf_high) -trapz (cpr_high) ) ; 

cd_low=s* (trapz (cpf_low) -trapz (cpr_low) ) ; 

edh=cd_high-cd; 

edl=cd-cd_low; 

ecd= (edh+edl) /2  ; 

cldoutput= [te (i) ’ , aoa (i) ’ , cl ’ , eel ’ , cd ’ , ecd ’ ] ; 
else 

%+l  degrees  AoA 
xt=xcpl (1) : S : xcpl (17) ; 
xb=xcpl (30) : s :xcpl (18) ; 
yc=-0.03490:s:0. 06406; 

%lift  coefficient  calculations 
%create  variables 
xct=xcpl ([1:17]); 
xcb=xcpl (  [18:3  0]  )  ; 

Cpt  =  Cp (  [1 : 17]  )  ; 
errt=eb ( [1:17]  )  ; 
cpt_high=cpt - errt ; 
cPt_l°w=cpt+errt ; 
cpb=cp (  [18 : 30]  )  ; 
errb=eb ( [18:30]  )  ; 
cpb_high=cpb+errb  ; 
cpb__low= cpb  -  e  r  rb  ; 

%spline  data 

f cpt=interpl (xct , cpt , xt , ’ linear ’ ) ; 

f cpt_high=interpl (xct , cpt_high, xt , ’linear*) ; 

f cpt_low=interpl (xct , cpt_low, xt , ’ linear ’ ) ; 

fcpb=interpl (xeb, cpb, xb, ’ linear ’ ) ; 

f cpb_high=interpl (xeb, cpb_high, xb, ’linear*) ; 

f cpb_low=interpl (xeb, cpb_low,xb, ’linear’) ; 

%calculate  CL  and  error  bars 

cl=s* ( -trapz (f cpt) + trapz (f cpb) ) ; 

cl_high=s* { -trapz (f cpt_high) +trapz (f cpb_high) ) 

cl_low=s* (-trapz (fcpt_low) +trapz (fcpb_low) ) ; 

elh=cl__high-cl  ; 

ell=cl-cl_low; 

ecl=  (elh+ell) / 2 ; 

%drag  coefficient  calculations 
%create  variables 


ycf =ycpl ([1:9,18:30]); 
ycr=ycpl ( [10 : 17] )  ; 
cpf =Cp (  [1:9,18:30]  )  ; 
errf =eb ( [1:9,18:30] ) ; 
cpf_high=cpf +errf ; 
cpf _low=cpf -errf ; 
cpr=cp ( [10 : 17] ) ; 
errr=eb ( [10 : 17] ) ; 
cpr_high=cpr-errr ; 
cpr_low=cpr+errr ; 

%spline  data 

f cpf  =  interpl (ycf ,cpf ,yc,  'spline  1 ) ; 
f cpf_high=interpl (ycf , cpf_high,yc, ' spline' ) ; 
f cpf_low=interpl (ycf , cpf_low, yc , ' spline ' ) ; 
f cpr=interpl (ycr, cpr,yc, ' spline ' ) ; 
f cpr_high=interpl (ycr, cpr_high,yc, 'spline') ; 
fcpr_low=interpl (ycr, cpr_low,yc, 'spline') ; 

%calculated  Cd  and  error  bars 
cd=s* (trapz (cpf) -trapz (cpr) ) ; 

cd_high=s* (trapz (cpf_high) -trapz (cpr_high) ) ; 

cd_low=s* (trapz (cpf_low) -trapz (cpr_low) ) ; 

edh=cd_high-cd; 

edl=cd-cd_low; 

ecd= (edh+edl) / 2 ; 

cldoutput= [te (i) ',aoa(i) ' , cl ' , eel ' , cd’ , ecd ' ] ; 

end 

if  aoa(i)==0 
if  te(i)==0 

cldout= [Re (i) ,cldoutput] ; 

fprintf (fidl, '%8i  %i  %i  %6.5f  %6.5f  %6.5f  %6.5f  \r\n ' , cldout) ; 
else 

cldout= [Re (i) (cldoutput] ; 

fprintf (fid2, '%8i  %i  %i  %6.5f  %6.5f  %6.5f  %6.5f  \r\n ', cldout) ; 

end 

elseif  aoa(i)==-l 

cldout = [Re (i) , cldoutput] ; 

fprintf (fid3, ' %8i  %i  %i  %6 . 5f  %6.5f  %6.5f  %6.5f  \r\n ', cldout) ; 
else 

cldout= [Re (i) , cldoutput] ; 

fprintf (fid4, ' %8i  %i  %i  %6 . 5f  %6.5f  %6.5f  %6.5f  \r\n ', cldout ) ; 

end 

%hold  on 

%plot (xt , -f cpt , xb, -f cpb) 

%clear  unused  variables 

clear  Rd  cd_high  cd_low  cljiigh  cl_low  cpb_high  cpb_low  cperrl  cperr2 
cpf  cpf_high 

clear  cpf_low  cpr  cpr_high  cpr__low  cpstdl  cpstd2  cptjiigh  cpt_low  edh 
edl  elh  ell 

clear  errb  errf  errr  errt  fcpbjiigh  fcpb_low  fcpf  fcpf_high  fcpfJLow 
f  cpr 

clear  fcpr_high  fcpr_low  fcpt_high  fcpt_low  he  noflowi  nstdl  nstd2 
ntap  ntdl  ntd2  nvalve 

clear  offsetl  offset2  stdl  std2  tau  tdl  td2  valve  ycf  ycr 
clear  fid  cpdata  m  cldoutput  cldout 

% NOTES : 

%Re ,  U,  aoa ,  inf  ile ,  inf  ilei  ,  nof  lov/,  nu,  out  file  ,  out  f  i  lei  ,  rho,  dl,  hd-  data 
% cp , cpl , cp2 , eb, ebl, eb2 , cl , eel , cd, ecd-  results 


% s ,  xcO , ycO , xcnl , ycnl , xcpl , ycpl , xc , yc , xct , xcb , cpt , cpb , f cpt , f cpb- 
useful  variable/data  vectors 

end  %main  for  loop 
f close (f idl) ; 
fclose (f id2) ; 
fclose (f id3) ; 
fclose (f id4) ; 
fclose (fid5) ; 
fclose (f id6) ; 


Appendix  G:  Data  Log  File 


This  appendix  contains  a  print  out  of  the  log  file  used  to  keep  track  of  all  data  associated 
with  each  individual  static  pressure  data  file,  such  as  flow  speed,  angle  of  attack,  trailing 
edge  configuration  and  water  temperature.  This  log  file  was  used  by  the  Matlab 
computer  code  to  aid  in  the  data  regression  used  to  analyze  the  data  presented  in  this 
document.  The  Matlab  computer  code  can  be  found  in  appendix  F. 
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